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ABSTRACT 



HIGGS BOSON PRODUCTION IN 
HADRON-HADRON COLLIDERS 

By 

Christopher James Glosser 

The Higgs boson is the last remaining component of the Standard Model of Particle 
Physics to be discovered. Within the context of the Standard Model, the Higgs boson 
is responsible for the breaking of the SUl(2) x U(1) gauge symmetry and provides a 
mechanism for the generation of the masses of the corresponding gauges bosons, the 
and the Z. In addition, this same mechanism provides masses for the leptons 
and quarks via Yukawa couplings. Therefore, the discovery and subsequent study of 
the Higgs boson and its properties is of the highest priority in particle physics today. 

The dominant production mechanism for Higgs bosons in hadron-hadron col- 
liders is via gluon-gluon fusion, in which gluons fuse via a virtual top quark to produce 
a Higgs. At the Large hadron CoUider, which has a center of mass of approximately 
14 TeV, one expects 4000-40000 of these particles to be produced via this process 
in one year of running. At leading order in Quantum Chromodynamics (QCD), the 
process gg ^ H produces a Higgs with no momentum component transverse to the 
beam axis. Higgs bosons of non-trivial transverse momentum may be produced if the 
the Higgs recoils against one or more emitted partons. 

In this work we seek to compute precisely the rate at which these Higgs par- 
ticles are produced as a function of the transverse momentum and the rapidity at 
the Large Hadron CoUider. This calculation is carried out to next leading order in 
the QCD strong coupling ag- In particular, this calculation includes the dominant 



gluon-initiated mechanism, which is responsible for the majority of the cross section. 
We work in the infinite top mass approximation, and argue that this is very precise 
in the domain of apphcabihty. 

We find an increase in the cross section by a factor of 1.9 - 2.0 in the interme- 
diate transverse momentum range, where this theory apphes. We also find that the 
renormalization scale dependence is much reduced in our calculation as compared to 
the leading order prediction. 
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Chapter 1 



Introduction: The Standard Model 
of Particle Physics 

The Standard Model (SM) of Particle Physics is the most stringently tested theoretical 
description of the world in which we live. It has succeeded in describing all subatomic 
phenomenology seen to date, to an incredible degree of accuracy However, even 
with all of the successes of the SM, a number of theoretical questions remain. Of these 
the most crucial is the mechanism in which the electroweak symmetry is broken. Over 
the past twenty years or so, a number of different methods have been proposed to 
generate this symmetry breaking. All models using a perturbative approach to this 
symmetry breaking require at least one additional scalar electroweak isospin doublet. 
The phenomenological manifestation of this scalar is generically dubbed a "Higgs 
boson." 

In order to understand why such a field is necessary, a pedagogical review of 
the Standard Model and its gauge theory is beneficial. From this discussion, we will 
come to understand the role of the Higgs field in maintaining gauge invariance, as 
well as understanding some of the theoretical problems that a lone Higgs doublet 
generates. 
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It is crucial to understand the experimental signatures of the Higgs with the 
next generation of high energy experiments. As the current phenomenological under- 
standing of the Higgs is plagued with large theoretical uncertainties, new calculations 
are required to improve the present predictions. With the discovery of the Higgs and 
a detailed theoretical and experimental description of its properties, we will gain our 
first insight into the question of electroweak symmetry breaking. 



1.1 Yang-Mills Gauge Theory 

The entire dynamics of the SM is the result of a particular manifestation of Yang 
Mills gauge theory 0, which was inspired by General Relativity. In Yang-Mills gauge 
theory, we note that measurable quantities are proportional to the absolute square of 
a complex field, and thus the general theory should be invariant under a general phase 
transformation. That is, we may demand that the fermion field ^(a;^) be invariant 
under the transformation 



Since from quantum mechanics, we know that all observables are computed from 
expectation values of some differential operator, it is quite trivial to see that these 
observables will remain unchanged under this phase rotation. This invariance is 
known global gauge symmetry. 

Now, we may require that the theory be invariant under a local gauge symme- 
try, in which the phase parameter is a function of the space time co-ordinates x^, by 
making the replacement 
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A typical fermion Lagrangian is of the form 

LF = ^>{x^){ip-m)^!{x^). (1.3) 

It is quite trivial to see that the Lagrangian is not invariant under such a trans- 
formation, since the derivative in the kinetic energy may act directly on the a{x^). 
However, if we replace the derivative with 

D^^d^ + igA^{x) (1.4) 

and demand that the field transform as 

A^{x) A^{x) + ^d^,{a{x^)), (1.5) 

wc see that the latter term cancels the anomalous term resulting from the differen- 
tiation of the phase. The phase parameter a{x) is, in this scalar function of 
the space-time variables. This is known as an abelian gauge symmetry. 

Now, we wish to know what will happen if we give the fermion field some 
gauge index, i, and demand that it transform under a local gauge symmetry 

^,{x) ^ [exp(ia«(a;,)T»)] ^,{x), (1.6) 

where we imagine that ^'i(x) transforms under some gauge group G. The matrices 
T" are the gauge group generators of some representation of the gauge group G. 
This is an example what is commonly called a non- abelian gauge transformation. 
Again, it is trivial to see that kinetic terms in the Lagrangian will be spoiled by this 
transformation. One then notes that the kinetic terms may be rescued, provided that 



3 



the conventional kinetic term id^ is replaced with the so-called covariant derivative, 



(1.7) 



This expression is completely analogous to equation (p,.4[). The gauge field v4^(x) is 



required to transform under the gauge transformation as 




(1.8) 



in order to cancel the term coming from the differentiation of the phase in the kinetic 
term. After making the replacement d ^ D, we find that the covariant fermion 
Lagrangian for a non-abelian gauge theory, 



is completely invariant under the transformation (|1.6| ). 

In order to have nontrivial dynamics, kinetic terms for the gauge fields must 
be introduced. These kinetic terms must satisfy the following properties: 

• They must be gauge invariant. 

• The resultant quantum theory must be renormalizable. ^ 

• There must be one term which is quadratic in d^. 

All of these constraints are satisfied by the tensor , where 



^ It is possible to construct models with "higher dimension operators" which are not renormaliz- 
able in the conventional sense. The divergent loop terms are then cut off at some arbitrary scale at 
which new physics comes into play. These theories are beyond the scope of the discussion, however, 
as they do not play a role in conventional Standard Model dynamics. 




(1.9) 



rpa pi 



ta 
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= T^{d,A:-d,A';^-gr"^AlAl}, (1.11) 
which allows us to write down a kinetic term for the gauge fields, 



Lgf = -|Tr[F^,Fn. (1.12) 



Note that a general gauge field is self-interacting because of the self-coupling term 
proportional to fate in the tensor F^^. f\ 

So far, everything is completely consistent, although we have been considering 
only the classical limit of the theory. If we wish to consider a Quantum Field Theory, 
in which the fields themselves consist of discrete numbers of excitations, we must be 
more careful in order to ensure that the unitarity of the S-matrix is not violated. There 
are two different procedures for quantizing a field theory, the canonical approach and 
the path integral ^ approach. After naive quantization of the theory, one finds that 
quantum corrections apparently spoil the unitarity of the S-matrix, for particular 
gauge choices. In order to compensate, one then adds ghost fields (scalar fields that 
obey fermion statistics) to cancel the non-unitary terms in the Lagrangian. 

While this may seem ad-hoc, the ghost terms may be derived using an explicit 
gauge fixing constraint in the path integral in what is known as Fadeev-Popov gauge 
fixing, in which the ghost terms are generated spontaneously when the gauge in the 
path integral is fixed 0]. Consider the path integral for a non-abelian field A^, 



/ VA exp [iS[A]] = JVA exp J d\ [{Gl^f) 



;i.l3) 



This path integral infinitely overcounts the field configurations, since it sums over 



2 



It should also be noted that a CP violating term ~\Tr[e^^pcyFp„F^^] is also consistent with 



the listed conditions. Again, this is beyond the scope of the discussion 

■^For an introduction to the path integral, a field theory text such as reference [|| is useful. 
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the gauge "copies" that represent the same physics. To remove this overcounting, we 
must estabhsh a gauge choice by constraining the field A through some functional 
relationship F[A] = 0. This is known as fixing the gauge. As an example, the familiar 
Lorentz gauge is fixed through the constraint F[A] = d/^A"''^ = 0. 

To fix our choice of gauge, we begin by inserting an expression for 1 into the 
path integral. 



Va6{F[Ao])Det 



'dF[A 
da 



where 



AT = A''^ + -D^.a" 



1 
9 



(1.14) 



;i.i5) 



We may freely change variables in the path integral from A A^. Since the inte- 
gration measure IDA and the action S[A\ are gauge invariant, everything in the path 
integral now is a function of A^ and therefore we may, without loss of generality, 
remove the label a. The resulting expression is 



VA exp [iS[A]] = ijva) J VA exp [iS[A]] 6{F[A])Det 



dF[A] 
da 



;i.i6) 



where the term / Va contributes an infinite overall constant to the normalization. 

The factor Det [dF[A]/da] can be evaluated by introducing the ghost field. It 
is instructive to work this out explicitly for the Lorentz gauge. Since the gauge field 
varies infinitesimally as (v4")" — (A^)'^ ^ -D^a"", the determinant may be rewritten; 



Det 



dF[A^] 
da 



Det 



-d^D„ 



Jvcjvc exp i J <fxc{-d^D^)c 



(1.17) 



:i.i8) 



where the c is an anti-commuting scalar field^, and therefore a ghost. We have thus 



^Anticommuting fields become Grassman variables in the path integral representation. 
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Boson Symbol Mass(GeV) Charge Color Isospin 

gluon g Adjoint No 

photon 7 No No 

W W± 80.41 ±e No Yes 

Z Z 91.187 No Yes 

Table 1.1: Gauge Bosons 

generated the Fadeev-Popov ghost term as a result of rigorously fixing the gauge. 
The c fields are indeed ghosts, and have no asymptotic states. They do, however, 
contribute as closed loops in diagrams. Moreover, the infinite gauge repetitions have 
been systematically eliminated in the path integral, and the gauge has been fixed. 

1.2 The Standard Model 

The Standard Model is, in essence, a Yang-Mills gauge theory, consisting of three 
"copies" of fermions, called fiavors, interacting with a set of gauge bosons dictated by 
the direct product of groups G = SUc(3) x SU/(2) x Uy(l). As was seen before, for 
each of these symmetries we may associate a dynamical gauge field. Therefore, there 
is a gluon field that interacts with the SU(3) "color" quantum number, A Weak 
Boson field WJ^^ that interacts with the isospin quantum number (/), and an abelian 
field that interacts with the Hypercharge {Y). 

1.2.1 Spin 1 Gauge Bosons 

Since the Electroweak symmetry is broken, we will observe admixtures of its fields, 
WJj^ and i?^, in nature. The 4 fields combine into two electric charge eigenstates, W±, 
a massive chargeless boson, the Z, and a massless chargeless boson, the photon 7. A 
list of the general properties of these particles may be found in Table |]l] 
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Fermion 


Symbol 


Mass(GeV) 


Charge 


Color 






nnn^i i 


—1 p 




muon 




.106 


-le 


Singlet 


tau 


T 


1.78 


-le 


Singlet 


ciGL^ li Uii ii" U. Li liiU 


7 / 


\ 1 tj /\ 1 u 


n 




muon neutrino 




< .00017 





Singlet 


tau neutrino 




< .0182 





Singlet 


up 


Ll 




-f 


1 "PI 1"^ 1 dT 

±1 ipiei 


charm 


c 


1.15 to 1.35 


|e 


Triplet 


top 


t 


175 




Triplet 


down 


d 


.003 to .009 


~x 


Triplet 


strange 


s 


.075 to .170 


Triplet 


bottom 


b 


4.0 tp 4.4 


Triplet 



Table 1.2: Standard Model Fermions 
1.2.2 Spin 1/2 Fermions 

There are three different generations of fermions. The first generation consists of 
the electron, the electron neutrino, the up quark and the down quark. The second 
consists of the muon, the muon neutrino, the charm quark, and the strange quark. 
The third generation consists of the tau, the tau neutrino, the top quark and the 
bottom quark. 

The particles of different generations may be grouped into subsets of particles 
which are indistinguishable, save their masses and Yukawa couplings (with the Higgs). 
These are known as families. For instance, an electron, muon and tau have exactly 
the same isospin and charge characteristics, and hence belong to the same family. 
Table |1.2| organizes the various known fermions by their family^. 



Since the left and right handed components couple differently to the elec- 
troweak sector of the theory, it is beneficial to list the separate components and their 



^The masses in Table 2 were obtained from reference The quark masses are the so-caUed 
"current-quark masses." 
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transformation properties in a separate table. The lepton doublets contain a left 
handed, -1 charge piece, which is "isospin down" (this refers to the eigenvalues of 
the (73 matrix) , and a charge zero component that has "isospin up" , to be identified 
with the e,/x,T (former) and their respective neutrinos (latter). We then arrange the 
components into isospin doublets. 



ye, J 



;i.i9) 



with identical expressions for the remaining two families. The Hypercharge (Y) may 
be computed directly from the charges and the isospin eigenvalue, using the for- 
mula Q — I + ^ (the charge here is in units of e). We find that the left-handed 
e,/x,T,fe,^'ju,and i/r have a hypercharge of -1, while the right-handed e,//,and r have a 
hypercharge of -2. 

We also arrange the left-handed quarks into an isospin doublet 



Qi 



(1.20) 



The left-handed quarks transform just as the left-handed leptons do under SU/(2), 
and have a hypercharge Y — ^. This, of course, implies that both components have 
nontrivial interactions with the electromagnetic field. The right-hand components 
of the up and down type components have hypercharges of | and — | respectively. 
Again, these may be computed directly from the charge formula, Q — I + ^. 

All quark species transform as triplets in the fundamental representation of the 
SUc(3) color force, while all leptons transform as singlets, and thus do not interact 
strongly. That is, each of the quark components may be written as a three component 
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J- K^L LILIU Ll 


a 


T 


J. 




Ll 


(-1,0) 


Doublet 


-1 


Singlet 


Ql 


(+2/3,-1/3) 


Doublet 


1/3 


Triplet 


Ur 


+2/3 


Singlet 


4/3 


Triplet 


dR 


-1/3 


Singlet 


-2/3 


Triplet 


Ir 


-1 


Singlet 


-2 


Singlet 



spmor, 



Table 1.3: Fermion Quantum Numbers 



Q 



;i.2i) 



where Q may be any of the quark fields, Qi, Ur or dr. The leptons do not couple to 
the color sector and hence are all singlets under SUc(3). The quantum numbers of 
the standard model particles are conveniently summarized in Table |1.3| . 



One may notice that, in dividing terms up into left and right-handed states, 
and by demanding that these states transform differently under the various gauge 
groups, that terms which mix the left and right handed terms are expressedly forbid- 
den. Since mass terms do exactly this, they are forbidden in (Unbroken) Electroweak 
gauge theory. Thus we must break the Electroweak gauge symmetry somehow in 
order to have a realistic model of the Standard Model Particles. 



1.3 Electroweak Symmetry Breaking 

The most glaring problem with the current particle physics phenomenology is the 
existence of massive particles. Electroweak gauge invariance forbids hard mass terms 
of the form /i^-B^, which is obviously not invariant under a gauge transformation 
B ^ B + dA. Also, the left and right-handed fermions, which transform differently 
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under SUi{2) and f/y(l), seem to preclude the existence of any mass terms for the 
fermions, since mass terms mix the right and left-handed components, once again 
spoiling gauge invariance. 

We wish to introduce these terms in a gauge invariant way. The first clue 
is that we must somehow include three additional degrees of freedom, one for each 
additional spin state of the three massive bosons. This seems to indicate that another 
field is necessary to compensate for the additional degrees of freedom. 

The general way to break the electroweak symmetry is through spontaneous 
symmetry breaking. To spontaneously break a symmetry, one introduces a field with 
a vacuum state which is not invariant under the symmetry in question. One then ex- 
pands the field around this vacuum. Since the kinetic energy terms in the Lagrangian 
of this field will, in general, contain covariant derivatives (which in turn contain cou- 
plings to the fields W,Z), expanding the field around its vacuum expectation value 
( VEV) , will introduce mass terms for the gauge bosons. 

If we also allow the field to interact with fermions via a fermion-antifermion- 
scalar interaction known as a Yukawa coupling, we find that we have also sponta- 
neously generated masses for the fermions. The mass of this field which we have 
introduced can be arbitrary (at least at tree level) in the SM. We therefore have a 
theoretical answer to all inconsistencies in the Standard Model. 

A scalar field which has these properties is generically dubbed a Higgs field. 
In the SM, the Higgs is a complex isospin doublet, so that it transforms under both 
the hypercharge and isospin rotations. Since this introduces four degrees of freedom, 
there is one left over degree of freedom after the Electroweak symmetry is broken. 
The phenomenological consequence of this is the existence of a massive, chargeless 
scalar called the Higgs Boson. 

One final caveat concerns the quark masses. The quark masses are not as well 
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defined as those of the other fermions, and so they must be treated with some care. 
Masses in a quantum field theory are basically due to the self coupling of the field. 
These, like other couphngs in the theory are affected by the renormahzation process. 
Since most masses in the SM are measurable parameters (that is, not predefined by 
the theory), particle masses are generally empirically defined as being the masses in 
the absence of other forces. The problem is of course that the quarks are confined 
inside hadrons and mesons, and thus are never "free particles" in the classic sense. 
For a quark, this limit never exists, so the idea of a "rest mass" is poorly defined. 
Quarks that have large Yukawa couplings to the Higgs may still have a rest mass 
defined, as long as it is larger than the strong coupling scale, h.qcD- 



1.3.1 The Higgs Boson 

To begin our discussion of the Higgs field, we write down the conventional Lagrangian 
for a standard model Higgs boson. It is 



2 ^'%2 



(1.22) 



where $ is complex scalar isospin doublet, and thus has hypercharge +1. The covari- 
ant derivative in this case is 



(1.23) 



where we include both the isospin field W^^ as well as the hypercharge field B^. Now, 
we parameterize the isoscalar doublet as 



V2 



;i.24) 
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or, writing this as a gauge rotation, 








v+H{x) 
\ V2 J 



(1.25) 



The dynamical field H{x) is the Higgs boson field. We see that for a particular 
gauge choice (commonly called the unitary gauge), we may eliminate three degrees 
of freedom. These degrees of freedom will show up as the necessary third degree of 
freedom for each of the gauge fields W^, S^, leading to the rather quaint analogy that 
the gauge fields have "eaten" three of the components of the Higgs doublet, leaving 
the Higgs boson field as an undiscovered dynamical field. 



As a consequence of choosing the field $ to have a nontrivial vacuum, numerous 
phenomenological characteristics of the SM emerge from the various couplings of the 
weak gauge fields with the Higgs field. Moreover, gauge invariant couplings with 
the quark and lepton fields generate mass terms for these particles as well, without 
breaking the overall gauge invariance of the theory. We begin by writing the most 
general parameterization for the model. The mass terms for the gauge bosons arise 
from the following term in the Lagrangian, 



1.3.2 The Masses of the W and Z Bosons 



{g2C7ijW;{x) + g^5,jB;ix)) (g^a^W^ix) + g^5j,B;{x)) 



(1.26) 



If we redefine the fields as 




(1.27) 



(1.28) 
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g^A^ = giW^ + g2B^, 



(1.29) 



where go — y 91 + 92-> '^^ fiiid that the following mass term arises: 

K (W^W;-) + \mI\Z,\\ (1.30) 

whereas the field ^4^ massless. This expression now contains explicit formulas for the 
gauge boson masses 

= gli^), (1.31) 
Ml = 9Uj)-{9l + 9l){jl (1-32) 

Hence, we identify the fields W^^ and as massive phenomenological gauge particles; 
whereas the field maintains a local U{1) symmetry and is identified with the 
photon field. It is straightforward to derive the expression for the electromagnetic 
charge e = by rewriting the covariant derivative in terms of the field A^. It is 
also customary to define the Weinberg angle, sin^^t^ = gi/ go- 



1.3.3 Fermion Masses and the CKM Matrix 

One may also postulate interactions between the Higgs field and the fermions of the 
standard model. We may write this Yukawa interaction with the fermions: 

+ (yrd'h^^Ql^ + iC.C.)), (1.33) 
14 



where n and m run over the fermion famihes. The field $^ is related to the complex 
conjugate of the Higgs field by 

$ = -i(Ty^*. (1.34) 

This expression is the most general term that we may write. Note that there is no 
requirement that the Yukawa matrix Unm be diagonal or real. It will, in general, mix 
electroweak eigenstates of the fermions. This has the implication that CP symmetry 
is violated in a general theory with a Higgs. 

First, let's consider the leptons. Since the neutrinos are considered to be 
massless, and have only a left hand component, we may freely redefine the fields 
such that the mass matrix is diagonal. The three diagonal coefficients of the mass 
matrices are then interpreted as the masses of the leptons. The diagonal Yukawa 
matrix elements are then the masses of the three leptons with the factor v/\/2 scaled 
out; that is, mi = yiv/\/2. Note that this also imphes that the Higgs boson couples 
to leptons with strength yi = \/2mi/v; that is, the coupling is directly proportional 
to the lepton mass. This is true of the Higgs couphng to all particles, and is of great 
importance in phenomenology, as we shall see in the next section. 

The up and down mass matrices, due to their additional complexity, may not 
be simultaneously diagonalized. It is conventional to choose a basis in which the up 
type quarks are diagonal, and the down type quarks are off-diagonal. One then finds 
that since the terms mix quark mass eigenstates, the W Boson violates "quark 
family number." 

Let's consider this situation in a little more detail. Since in the Standard 
Model the only coupling which actually mixes isospin-up and isospin-down fermions 
is the coupling to W^, we find family mixing terms in the couplings between the 
fermions and the W boson terms. That is, the W couples to "weak eigenstates" 
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which differ from mass eigenstates by the rotation matrix, V: 

I \ 



Weak 



\ 



Vud 




Vub 








Vtd 


Vts 





I 



:i.35) 



Mass 



This is the so called Cabbibo-Kobayashi-Maskawa (CKM) matrix ||5|. This matrix 
has the standard parameterization 



C12C13 S12C13 siacis 

— ■512C23 — Ci2C23Si3ei3 C12C23 — Si2S23Si3ei3 S23C13 
S12S23 - Ci2C23Ci3ei3 -C12S23 - Si2C23Si3ei3 C23C13 



1.36) 



where c,- 



cos [Pij),Sij 



sin(6'jj), and Cij = exp{—i6ij). In the Standard Model, 



these angles are not constrained by any symmetry, and must be determined experi- 
mentally. It may not be possible to make these same arguments in extended models, 
and one may consider searching for deviations in the CKM matrix as an indication 
of physics beyond the SM. See Q for an excellent discussion of these indicators. 

If the neutrino sector of the theory contains a Dirac mass term, we must forego 
the previous analysis, as one cannot simultaneously diagonalize the electron and neu- 
trino mass matrices unless at least one of the two is completely massless. Recent 
measurements of neutrino oscillations [0 do seem to suggest that the neutrino species 
have a mass difference (and therefore, masses), which of course must be accounted for 
theoretically. One may postulate the existence of an additional right-handed compo- 
nent, which will in general have its own coupling to the Higgs. After diagonalizing 
the charged leptons, the neutrino components are left with off-diagonal terms. These 
terms, while extremely important to astrophysics, do not contribute much to con- 
ventional high energy physics collider experiments, where the approximation that the 
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neutrino is massless is certainly viable. 



17 



Chapter 2 

The Higgs Boson: 
Phenomenological Overview 

The Higgs boson is currently the only missing component of the Standard Model. 
Precise knowledge of Higgs physics will most likely give the first glimpses of the un- 
derlying structure responsible for electroweak symmetry breaking. Since the couplings 
of the Higgs to the various particles in the Standard Model are already determined, 
any deviation from the rates computed using these couplings could be an indication 
of new physics. 

However, in order to examine these rates, we must first find a Higgs and mea- 
sure its mass. The Higgs mass is not fixed by any theoretical argument in the Stan- 
dard Model. In fact, a tree- level analysis quickly reveals that there are no constraints 
whatsoever on the Higgs mass. 

At the 1-loop and higher orders, however, constraints do arise on the Higgs 
mass. A renormalization group analysis reveals that the Higgs mass cannot be too 
large or too small relative to the Weak scale. The former results in triviality, while 
the latter in vacuum instability. The Higgs mass also appears in 1-loop corrections 
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to ee ^ Z ^ ff, and so by precisely measuring the cross section to this process, 
we may constrain the Higgs mass indirectly, as was done for the top quark before its 
discovery in 1995 g, |. 



We begin our analysis of Higgs phenomenology by considering some theoretical issues 
concerning the mass of the Higgs. As mentioned before, the classical theory of the 
Higgs and its couplings to the other Standard Model particles does not in any way 
restrict the Higgs mass. However, radiative corrections cause the VEV and the quartic 
coupling to run, and can cause problems if the Higgs mass is too far removed from 
the other electroweak parameters. 

One-loop corrections to the Higgs effective potential lead to vacuum stability 
constraints. These result from the fact that the asymmetric electroweak-violating 
vacuum could decay into a symmetric one, therefore destroying the universe as we 
know it. Typically, these analyses put lower bounds on the Higgs mass which are far 
less then the current mass bound from LEP. 

Triviality arises from the fact that the quartic coupling A in 0^ theory exhibits 
a Landau pole [0, and therefore the theory becomes strongly coupled at some high 
energy scale. To see this, consider the expression for A running from some energy 
scale /io to /i. 



Since Aq is related to the Higgs mass through the VEV (i.e. rrijj = Af^), the re- 
quirement that the coupling remain finite below some relevant scale puts an upper 
constraint on the Higgs mass. If the theory is to remain perturbative up to the Planck 
Scale, then the Higgs mass cannot be more than about 1 TeV. 



2.1 Vacuum Stability and Triviality 




(2.1) 
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Since this behavior is essentially nonperturbative, such calculations are best 
done on the lattice. [] A recent lattice study ||Tl| has placed an upper bound on the 
Higgs mass of about 460 GeV, if no new physics arises before then 0. 

2.2 Electroweak Precision Data and the Rho Pa- 
rameter 

The masses, decay widths and cross sections for W and Z also receive corrections 
from virtual Higgs loops. The result is that the p parameter, 



which is 1 at tree level, picks up the following corrections from virtual top and Higgs 
loops: 



All of the fermion loop corrections to p vary quadratically with their masses , so 
it is easy to see that the top gives the dominant contribution. Therefore, a precise 
knowledge of 9w ^w, f^z, and mt can be used to constrain the Higgs mass. 

The major problem with indirectly measuring the Higgs mass is that the cor- 
rections to the p parameter vary only logarithmically with Mh, rather than quadrat- 
ically, as do the top corrections. It is hence very difficult to have enough statistics 
and eliminate enough systematic uncertainty to put a precise constraint on the mass. 

As of this writing, the upper limit from the electroweak precision data is 224.8 

^Lattice Gauge Theory is an attempt to formulate quantum field theories numerically by placing 
them on a discrete space-time lattice, and evaluating the path integral numerically. 

^Reference is an excellent introduction to these theoretical issues, including the theoretical 
issues surrounding the Higgs phenomenon in extended models. 
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Figure 2.1: Higgs Width vs Mass 

GeV. The lower bound on the Higgs mass is most strongly constrained by the direct 
search at LEP|12[- They find a lower bound of 114.1 GeV at 95% confidence level. 



2.3 Phenomenology at Colliders 

In order to search for the Higgs in a collider experiment, we must have a precise 
understanding of how it is produced, and how it decays. This is the main topic of 
this section. 

2.3.1 Width and Branching ratios 

A graph of the Higgs width as a function of mass is provided in Figure |2]^ 0. As 

can be plainly seen, the width increases dramatically around 160 GeV. This dramatic 

^This graph and the branching ratio graph (figure |2.2| ) were generated using the program HDE- 
CAYlH 
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increase in the width is due to the fact that the Higgs is now massive enough to decay 
into weak bosons, with which it has a comparatively strong couphng. Higgs particles 
that can decay into weak bosons have a collider phenomenology dramatically different 
than those which cannot decay to weak bosons. The former are referred to in this 
work as heavy Higgs, and the latter intermediate mass Higgs. 

Looking at the branching ratios for a Higgs in the intermediate mass range 
(Figure p^) , we see why its detection is so problematic. As noted in the last section, 
the Higgs couples to particles with a coupling strength that is proportional to the 
particle mass. Therefore, it prefers to decay to the heaviest particles for which it is 
kinematically allowed. For the intermediate mass Higgs, it decays most of the time 
into bottom quarks, which are swamped by QCD background events. If one wishes 
to look for the Higgs through the bb channel, one must be able to efficiently tag 
bottom quarks, via vertex detection. It is not clear that bottom tagging will have the 
resolution necessary to see Higgs events. 

About 18% of the time, the Higgs will decay into taus, charms or gluons. The 
decay to taus is also a very promising mode, since taus usually decay into electroweak 
particles. Observations of Higgs decays to charm quarks is basically out of the ques- 
tion, since the ability of detectors to tag charm quarks is extremely unreliable. Gluons 
are out of the question as well, since they may not be tagged at all. 

The mode H —>■ 77, which arises through loop effects, is quite rare, with a 
branching ratio of 10~^. This, in particular, is quite unfortunate, since this mode 



offers the best chance for detection in a hadron-hadron collider [|14[ [1^. This process 
has been called the "golden mode", since the energy resolution of these events in a 
particle detector is relatively fine and its background is fairly small. One can hope to 
measure the rate to 7's and look for the characteristic enhancement in the transverse 
momentum spectrum, provided that one has a precise prediction of the spectrum of 
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Figure 2.2: Higgs Branching Ratios 

background events. 

The situation is better for a heavy Higgs, which decays predominately to 1^'s, 
Z's or top quarks. While the production rate for a Higgs in this mass range is just 
too small to measure at the Tevatron, these rates can presumably be measured at the 
LHC, provided that the Higgs mass is within this range. 

2.3.2 Cross Sections 

Overall, the best way to search for a Higgs is with a lepton collider, such as LEP, since 
one has a precise knowledge of the initial state energies of the interacting particles. 
At these machines, one begins to expect Higgs particles around an energy of £'com = 
Mz + Mh, with the Higgs radiating from the final state Z boson. The decay into 
jets or leptons is relatively clean, allowing the invariant mass of such events to be 
reconstructed. If the Higgs mass were within the accessible range of LEP, then it 
most surely would have been discovered. This allows a lower bound to be placed 
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on the Higgs mass, which is currently 114.1 GeV. Recently, a recent analysis by the 
ALEPH collaboration revealed an excess in events at LEP at around 206 GeV [[12 1. 



These events, while still preliminary, are very suggestive of a Higgs with Mh ~ 115 
GeV. 

Producing a Higgs boson at a hadron-hadron collider is generally not difficult. 
In fact, even at the Tevatron (FNAL), there are most likely dozens of Higgs produced 
each year that the laboratory takes data (depending on the Higgs mass). The real 
problem is detection; if the mass is less than 135 GeV, the Higgs decays predominately 
into hb pairs. Since these particles generally show up in the detector as jets, this 
process is swamped by QCD background events. Therfore, for Higgs bosons below 
135 GeV, the most likely channel for discovery at the Tevatron is the associated 
production channels pp — HZ and pp — > HW . For these channels, tagging of b- 
quarks is crucial. For itlh > 135 GeV, the WW* decay mode (where W* represents 
an off-shell W boson) becomes available. The actual reach at the Tevatron at Run H 
depends on the total integrated luminosity, of course. A recent study |jT6| found that 



for 20 fb \ the Tevatron can exclude a Higgs at 95% confidence level up to 180 GeV, 
and could discover it at the 3a level if the mass is below 125 GeV. 

If the Higgs boson is not discovered at the Tevatron, then the Large Hadron 
Collider (LHC) at CERN, which is due to start runs around 2006, will find the Higgs 
particle. The main advantage of the LHC is that there will simply be more events, 
due to the extremely large luminosity of this machine, as well as the enhanced cross 
section. This allows for a much better chance for detection, especially in the if — > 77 
mode. We will therefore focus our discussion on this device, where the discovery of 
the Higgs will most likely happen. 

There are a number of ways that a Higgs boson can be produced at a hadron 
collider. One notable mechanism is the production rate through the gluon fusion 
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Figure 2.3: Higgs production rates at the LHC 

process, gg H, which dominates the overall production rate. Other ways to pro- 
duce a Higgs generally fall into the category of radiation from electroweak bosons, 
production through electroweak boson fusion, and radiation from heavy quarks. The 
expected values of these cross sections at the LHC are given in Figure ^]3| Q 

Clearly, the process gg ^ H dominates the total production rate at the LHC. 
Therefore, we will concentrate our discussion on this particular production mecha- 
nism. This mechanism is not only the most dominant production mechanism, but 
is the least understood in a quantitative sense. In particular, we will show that the 
QCD corrections to the overall cross section as well as the differential cross section 
are sizeable. 



"^We thank Michael Spira for the generous use of this plot, as well as the program HDECAY 
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Chapter 3 

Perturbation Theory 



If one wishes to search for a Higgs particle in a high energy coUision, one needs to 
have an accurate prediction of its behavior in such an event. For instance, we would 
like to know in what kinematic regions we should should search for decay products. 
In order to do so, we are going to have to calculate observables using the theoretical 
apparatus of Quantum Field Theory. 

Integrations over momentum in such theories are plagued with divergences of 
several types. We will therefore have to come up with a consistent way to render 
these integrals finite, so that the divergences may be dealt with in a rigorous manner. 

Once the divergences are regularized, we need a systematic procedure for elim- 
inating them from our predictions. Fortunately, in a well defined theory, these diver- 
gences will either cancel amongst themselves or we will be able to factorize them into 
universal parameters of the theory in a process known as renormalization. These 
procedures, if applied rigorously and correctly, will allow us to compute precise, 
divergence-free expressions for observables in particle accelerators. 
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3.1 Divergences and Regularization Schemes 



3.1.1 Divergences 

When one calculates an amplitude in QCD or any other field theory, one typically gets 
expressions which are poorly behaved in the large energy region of phase space. These 
are known as ultraviolet (UV) divergences, since they correspond to a divergence in 
the number of high energy (frequency) fiuctuations. An example of an integral which 
contains these types of divergences is 



Clearly this expression diverges in the limit that the regularization parameter A goes 
to infinity. 

In addition, if the mass of the particles is small, two additional types of diver- 
gences appear, which are collectively called infrared (IR) singularities. First, there is 
a proliferation in the number of particles generated at low energy. This is known as 
a soft singularity, called such since it causes little recoil when emitted from a more 
energetic particle. In a renormalizable quantum theory such as QCD, these infrared 
singularities will cancel in the calculation of any properly defined "infrared safe" ob- 
servable quantity. We will return to the question of what it means to be an "infrared 
safe" observable later. 

The other type of infrared singularity that one may have is called a collinear or 
mass singularity, which occurs, for instance, in QCD. A collinear singularity occurs 
when the momenta of two massless particles line up. To be more precise, as the 
invariant mass of the two distinct particles becomes zero, the resultant jets that are 
produced when the particles hadronize become indistinguishable from one another. 
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Therefore, one does not really know whether a jet of hadrons came from one particle 
or two or more that were collinear. We will discuss how to handle such divergences 
in the section on perturbative QCD. 

3.1.2 Historical Regularization Schemes 

In this subsection, we shall discuss two regularization schemes that were used histori- 
cally to quantify the ultraviolet and infrared divergences . We will see that they both 
lead to consistency problems. 

In the fictitious mass scheme, one introduces a fictitious mass term into the 
integration to regularize the soft and collinear divergences. The problem with this is 
that since mass terms violate gauge invariance, gauge anomalies will in general arise 
in calculations, and their removal is quite subtle. For this reason, this fictitious mass 
scheme is rarely used in modern calculations. 

The UV divergences were historically handled using the Pauli- Villars regular- 
ization scheme (PV). PV regularizes the UV divergences through the introduction of 
a massive scalar ghost field, which converts the problem terms into logs of the PV 
mass term Apy. The ghost fields must be introduced in a gauge- invariant manner, 
which can become quite complicated for all but the simplest theories. 

3.1.3 Dimensional Regularization 

Almost all current calculations use dimensional regularization (DR) to regularize the 
divergences. Using dimensional regularization has a number of advantages. Since no 
extraneous fields are introduced, unitarity is manifest, and the gauge invariance of 
the theory is preserved. 

Let us explore the basic ideas of dimensional regularization. In a divergent 
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field theory, we expect to have terms arise in calculating observables that are poorly 



behaved in the large (UV) or small (IR) energy limit. Equation contains examples 
of both of these types of divergences First, it diverges as A —>■ oo. This is the 
aforementioned ultraviolet divergence. Second, the integral also contains an infrared 
divergence, as m ^ 0. 

The main idea of dimensional regularization is that divergent integrals of this 
type may be rendered finite by analytically continuing the number of space-time 
dimensions to 4 — 2e dimensions. The ci(p^) terms which come from the radial piece 
of the phase space in 4 dimensions, are replaced with the analogous piece in 4 — 2e 
dimensions: 

d^p = ^d{p^)dT ^ d'-^'p = (3.2) 

Where dV^ is the 4 — 2e dimensional angular measure. The net result of our efforts 
is that the divergent p"^ integral, which arises from evaluating the phase space in 4 
dimensions, picks up an additional factor p"^"^: 



i 



^P'^l", (3.3) 



— \ A*^ / ' 

which is well defined for e > 0. The n term, for now, maintains the dimensionality of 
the original expression. 

We may also regularize the IR divergence in a similar manner. Eliminating 
the mass term in equation ( p.3|) , and noting that we must carefully separate the IR 
divergent and the UV divergent regions of phase space, we obtain 
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where we have defined euv > and ejR < to represent the uhraviolet and infrared 
poles, respectively. In QCD, we also find that poles arise from the angular integration 
piece, where the internal partons become coUinear. These poles in e also contribute 
infrared poles to the matrix element. 

In dimensional regularization, then, we find that we have traded a poorly 
behaved log for a pole of the form, — i. These poles must either cancel amongst 
themselves (which is the case for the infrared poles), or they must cancel against 
renormalization counterterms (as is the case for ultraviolet or collinear poles), which 
we will discuss shortly. If everything works as it should, we should be able to analyt- 
ically continue the complete expression, including counterterms, from e = ±00 
without hitting any poles. That is, all poles of the form (— e — n)^^, where n is a 
integer must cancel Q We may therefore set e — with good conscience. 

Now that we have examined the basic ideas of dimensional regularization, 
let us examine the finer details. Immediately, we see an ambiguity in this scheme. 
It is not immediately obvious how to handle quantities that are only well-defined 
in integer dimensions, like the number of spins, nor is it obvious what to do with 
antisymmetric tensors like 7^ (which depends crucially on the integer nature of Ndim)- 
The most common prescription is to continue the number of spins of bosons by setting 
Nspin = 2(1 — e), and by demanding that: 

• 7^ anticommutes with 7°, 7^, 7^, 7^, and 

• commutes with all other 7's. 

This convention was developed by t'Hooft and Veltman [T^, and is the standard way 



of handling this problem. Gauge invariance is manifest in this scheme, and there 



is the degree of divergence of the theory, that is, a quadratically divergent theory will 
typically have terms like (e — All of the interactions of the Standard Model save that of the 

Higgs are logarithmically divergent, so all of these superficially higher order divergences must cancel 
amongst themselves. 
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are precisely two chirality states (defined in terms of projection operators) for the 
fermions. Thus, this scheme avoids any anomalous degrees of freedom which could 
arise with other "generalizations" of 7^. By setting the number of gluon hehcities 
equal to Ngpin 2(1 — e), we avoid having to separate out the g^y tensors that arise 
from the internal spin sums. This way of handling the helicities, where the gauge 
boson spins are manifestly continued to 2 — 2e and the number of fermion helicities 
is fixed to be 2, is known as Conventional Dimensional Regularization. 



3.2 Renormalization 



Having regularized the divergent pieces, we now seek to eliminate them from our 
predictions in a manner which is consistent with our theory. We will take advantage 
of the freedom that our theory provides to absorb these divergent terms into the 
definitions of the wave function normalization, the definition of the charge, and the 
definition of the mass. That is, we have the freedom to fix our normalizations to 
whatever we wish. We may, for example, make the replacements 

e — > Zf.e, (3-5) 



These reparameterizations will allow us to absorb the ultraviolet divergences into 
these definitions. Provided that we do this consistently, we arrive at a well-defined, 
finite prediction for our observable. The gauge invariance of the theory will in general 
impose relations between these normalization constants, but in general, we still have 
significant room to absorb the ultraviolet divergences. 
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3.2.1 Modified Minimal Subtraction 



Even with the gauge restrictions, we find that we have significant freedom to repa- 
rameterize the theory. We must therefore set down conventions for the inclusion of 
the UV counterterms. These terms, which arise because of a redefinition of the nor- 
mahzation constants in the field theory, are somewhat ambiguous. This ambiguity 
arises since renormalization group analysis only defines how one modifies the charge 
to 0[e~^]. We wish to remove the poles in e, but we also have the freedom to add any 
finite piece to the renormalized charge that we wish. Consider 



where the coefficient (3o is the first renormalization group coefficient. The details of 
this coefficient depends on the symmetry of the group, and the number of fermions 
interacting with the gauge bosons. 0[e^] represents terms that are unspecified by any 
gauge or other restrictions. 

In modified minimal subtraction (MS), we keep finite terms which tend to arise 
from the finite phase space integration in the dimensional regularization procedure 
outlined above. Instead of subtracting just the e pole, we subtract off the piece 



in MS. Again, this is a just a convention which is used to remove the extraneous 
universal pieces that tend to complicate calculations in dimensional regularization. 
There are a number of other ways of expressing the MS counterterm, which are 
technically different from this one at 0[e]. However, since these 0[e] terms will 
eventually be set to zero in a next-to-leading order calculation, the same result is 
obtained. 




(3.6) 




(3.7) 
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In principle, the subtraction scheme should not matter when one is calculating 
an observable. In practice, however, different subtraction schemes can cause signifi- 
cant variations in the final result at a fixed order in perturbation theory. One must 
also consistently use the same subtraction scheme to calculate the parton distribution 



functions from the raw deep inelastic scattering (DIS) data [T^, as well as to calculate 
the UV counterterm from the charge renormalization. If one wishes to do a calcu- 
lation in a different scheme, one must either know how to convert the result in that 
scheme back into MS, or one must be able to show that the scheme produces the same 
results at that order of perturbation theory. Since it is fairly involved to do either 
of these things, most calculations are done using the modified minimal subtraction 
scheme. 



3.3 Quantum Chromo dynamics 

Quantum Chromodynamics is the theory of the strong interaction at high energies. 
This model introduced dynamics to the quarks of the Gell-Mann SU(3)flavor model, 
and proposed a carrier of the strong nuclear force, the gluon. We shall begin by 
discussing the Naive Parton Model, which explains the behavior of the hadrons in a 
qualitative sense, and then move on to perturbative QCD, which allows one to compute 
much more precise expressions for collider observable quantities. 

3.3.1 The Naive Parton Model 

We shall begin the discussion of Quantum Chromodynamics (QCD) by introducing 
the parton model, due to Bjorken and Feynman. The parton model was introduced 
in an effort to explain the deep inelastic scattering data, which seemed to indicate 
that the proton had a definite substructure. As originally proposed, the parton model 
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suggested that hadrons consisted of massless point particles which could then interact 
electroweakly with the scattering electrons to produce jets (semi-coUinear collections 
of hadrons). The model originally had no other interactions built in to explain any 
of the strong coupling effects seen. This model is known as the naive parton model, 
and ignores any perturbative QCD effects. Indeed, this model did seem to describe 
the data in a qualitative fashion, and it was later deduced from the DIS data that the 
components of the proton that were interacting with the electron were spin 1/2 ob- 
jects. This is precisely what one would expect it these "partons" were to be identified 
with the quarks of the Gell-Mann SU(3)flavor model. 

The naive parton model assumes the high energy limit, in which both quarks 
and hadrons are massless. To calculate a cross section in DIS, we convolve a partonic 
cross section with a form factor fff^qiC)- In other words, the hadronic cross section 
is 



where we have introduced fn-^qix), the process-independent parton distribution func- 
tions (PDF). These empirically measured functions represent the probabihty of the 
interacting parton having a fraction ^ of the parent hadron's momentum. In the 
limit that both particles are massless, the hadron momentum is given in terms of the 
proton momentum as 



The hadron-hadron cross section is defined in a similar fashion, where the 
structure functions are deemed to be universal, that is, process independent. We 
may then measure the parton distribution functions through DIS and apply them to 
any process we wish. The hadron-hadron cross section is then given by 




(3.8) 



(3.9) 




(3.10) 
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In equation (|3.8| ), q may be a quark or an anti-quark, and in equation (|3.10|) , a and 
h may be any species of parton. These expressions must be summed over the parton 
species, g, a, and 6, since partons may not be directly observed. 



3.3.2 Perturbative Quantum Chromodynamics: Improved Par- 
ton Dynamics 

Ideally, one would like to extend the qualitative successes of the naive parton model 
with a dynamic theory. Such a theory would have the same qualitative behavior of the 
naive parton model, but with enough computational power to make precise analytic 
predictions of observable quantities. 

We find that, if we allow the quarks of the Gell-Mann model to interact dy- 
namically via an unbroken SU(3) Yang-Mills theory, we succeed in extending the 
parton model. This theory contains generalizations of the hadron structure functions 
mentioned in the previous subsection, and allows improvements of the predictions of 
that model by the dynamic interaction with the SU(3) gluon field. 

This model based on the SU(3) gauge theory is called Quantum Chromody- 
namics (QCD). QCD is an (unbroken) SU(3) gauge theory in which the quarks of 
Gell-Mann's flavor model interact with spin-1 bosons known as gluons. There are 
eight of these, interacting with the quarks via a three-dimensional quantum number 
known as color. One may show that such a theory satisfies asymptotic freedom, in 
which the coupling becomes very small at high energy scales. 

This theory is also interesting, since it gives an explanation as to why the 
quarks are never seen in nature. The coupling as of the SU(3) gauge theory actually 
diverges for a small energy scale, indicating that any perturbation theory is poorly 
defined in this limit. In quantum mechanics, small energies are equivalent to large 
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distances. This implies that any attempt to pull quarks away from one another results 
in a steady increase in energy, until the available energy is enough to create a quark 
anti-quark pair out of the vacuum. In order to examine such behavior, one must be 
able to calculate exact expressions, without resorting to perturbation theory. 

If the interaction is over a very small distance (large energy scale) , one may use 
perturbation theory, since as is small over these distances. Such an expansion, and 
the requisite renormalization of the PDF's and the coupling constant, are known as 
perturbative QCD. For very "hard" scattering processes, where a significant amount 
of energy is exchanged among the hadrons, we may calculate corrections to the naive 
parton model. This consists of calculating a parton-level cross section in standard 
perturbation theory, then inserting the results into the convolution equations above. 

QCD, like all gauge theories, is poorly behaved in the ultraviolet region. Wave 
function renormalization removes the ultraviolet singularities that occur. In addition, 
two other types of singularities occur, collinear (or mass) singularities, and soft (or 
infrared) singularities. The soft singularities must cancel when all of the pieces of the 
cross section are added together, but the collinear singularities do not cancel against 
any other pieces of the calculation. However, the collinear poles are universal (process 
independent), and may be absorbed into the (universal) structure functions. That is, 

fh-.p{z) ^ I dx/xfh-^p'{x/z,iiF)Kp,^p{z), (3.11) 

t/ 

where the convolution kernel Kp^p/{z) is 

W.) ^.(1-.,.^- (i^) (l=^)'S=£M + fi„,e. (3.12) 

The convolution kernel may be thought of as being analogous to the normalization 
constants Z in the ultraviolet renormalization procedure. The splitting kernels P{z) 
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in equation (|3.12|) are 



F„W = + (3.16) 



In the preceding relations, we have introduced the "+" prescription, which is defined 
through the relation 

f = f (3,17) 

JO (1 — z)+ Jo (1 — z) 

for any arbitrary function g(z). 

In equation ( |3.12|) we have introduced another dimensional byproduct of the 



renormalization process, the factorization scale fip- We have also made explicit the 
separate scales introduced in ultraviolet and infrared regularization processes; the 
ultraviolet renormalization scale will be referred to as /i/j. In practice, one generally 
takes these scales to be the same, though there is no theoretical argument that they 
should be taken to be the same. When we have set these scales equal to each other, 
we will refer to the resulting parameter as fi, to emphasize it's generic nature. 

The physical interpretation of the factorization scale is somewhat different 
than that of the renormalization scale. Indeed, the factorization scale represents the 
cutoff scale for so-called "hard processes." In other words, it quantifies exactly how 
much of the soft QCD radiation is to be included in the partonic matrix element, and 
how much is to be absorbed into the parton distributions. The renormalization scale 
represents the cutoff scale at both the large energy end, where we expect new physics 
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to come into play, and the low end, where we absorb the remaining IR singularity 
into the quantum mechanical excitations. 

Now, in order to cancel the soft poles that occur in the matrix element, we 
must carefully define our observables. In general, an observable O may be a function 
of the number of particles, as well as the kinematics of the system. In order for an 
observable to be infrared safe, we guarantee that (for N and N + 1 parton systems): 

Oiv+i(£^i^0) = Ojv. (3.18) 

That is, 0{N + 1) is infrared safe if, as the energy of one of the partons vanishes, 
0{N + 1) becomes identical to the analogous function with one less parton, 0{N). 
Similarly, an observable is collinear safe if for two particles with momenta pi, p2, we 
have that 

On+i{Si2 ^ 0) = On, (3.19) 

where (5*12 = {pi +^2)^)- If O satisfies these two conditions, then the cross-section is 
guaranteed to be finite to NLO in perturbative QCD. 

This process may be generalized. We may calculate cross sections and other 
observables to any order in perturbation theory that we wish. These calculations tend 
to be quite comphcated, as the additional complexity of the SU(3) gauge group tends 
to generate large numbers of contributing graphs. This is especially true in processes 
that involve large numbers of gluons, since gluons can couple not only to fermions, 
but to other gluons as well. 
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3.4 Helicity Amplitudes and Color Ordering 



In recent years, many techniques have been developed to simphfy the computation of 
multi-parton amphtudes in QCD By using techniques such as hehcity am- 

phtudes and gauge-invariant color orderings, one may simplify a calculation greatly. 
This technique significantly reduces the number of diagrams that one has to evaluate 
in computing, for instance, gluon amplitudes. 

In order to avoid calculating things redundantly, one may make use of some 
of the symmetry properties of QCD. In the approximation that all of the partons are 
massless, the spinor states of the partons become true helicity eigenstates. Since some 
of these states are related to other states by charge conjugation, Bose symmetry, or 
parity, one avoids calculating redundant helicity states. Also, each non-cyclic color 
projection is gauge invariant from other non-cyclic projections, so we may reduce the 
number of terms by calculating a single cyclically related color projection. 

Another simplification lies in the fact that these projections of color and he- 
licity are gauge invariant, so that one may choose a different gauge to evaluate the 
partonic matrix element for each different helicity state. It is generally convenient to 
express polarization vectors in terms of massless Weyl spinors [^, |22| . The polariza- 
tion vectors for a gluon of momentum p may be written 

, , (±p\Y\q±) 



where 



\p±) = ^^u{p) (3.21) 



is a massless spinor, and q is some arbitrary massless reference momentum such that 
q ■ p 7^ 0. Since a change in q is equivalent to adding a piece proportional to p, we 
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have basically absolute freedom to choose q. The terms proportional to p will then 
cancel in gauge invariant amplitudes via the Ward identities. 

Using these expressions, the number of terms one needs to calculate a QCD 
observable is reduced considerably. However, ambiguities arise at NLO, where we 
would like to analytically continue the number of dimensions. The number of spins 
states in a conventional MS calculation is generally taken to be 2(1 — e) which is 
quite clearly an abstraction. If we consider actual helicity states, the number of 
external particles must be two exactly. Internal (loop) particles may have either, 
but are generally taken to have 2(1 — e) helicity states. Therefore, when calculating 
expressions using helicity states, one must know how to alter the resultant expressions 
in order to obtain the conventional MS expression. We will encounter this issue later, 
during our calculation of the Higgs differential cross section at NLO in QCD. 
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Chapter 4 



The Higgs Boson Cross Section and 
Transverse Momentum Spectrum: 
Leading Order Calculation 



One of the most important quantities that will be measured at the LHC is the pro- 



duction rate of the Higgs |Tj, |23|, ^ . Therefore, a precise theoretical understanding 



of the rate of Higgs production is critical to any attempts to search for the particle. 



Since the cross section is dominated by the gluon- fusion process this process 
must be calculated as precisely as possible, in order to precisely determine how many 
of these particles will be produced in a year of running. 

Aside from having a good idea of the overall production rate of the Higgs, the 
transverse momentum spectrum and the rapidity spectrum of the Higgs boson is of 
great interest to the LHC project. Since the detectors themselves must ideally be 
designed around the kinematic regions where one expects to see the Higgs, having a 
reliable prediction for these differential cross sections is of particular importance. 

Unfortunately, QCD observables are notoriously unreliable at leading order. 
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and these must generally be calculated to NLO (at least) in order to have a reli- 
able prediction for their behavior. Moreover, leading order QCD calculations tend 
to exhibit a strong dependence on the factorization/renormalization scale, fi, which 
reflects this theoretical uncertainty. 

The focus of this work is to calculate a reliable numerical estimate of the trans- 
verse momentum/rapidity spectrum. In order to improve the existing calculation, we 
will add to it radiative corrections which will hopefully have less scale dependence, 
and therefore less theoretical uncertainty. Since the procedure for calculating this 
differential cross section at NLO includes calculating the LO spectrum, it is instruc- 
tive to begin with the much simpler LO calculation. We will then show in the next 
chapter how to modify this expression in order to obtain the differential cross section 
at NLO. 

We will assume in this calculation that the dominant contributions to the 
transverse momentum spectrum come from the aforementioned gluon-fusion process. 
We will calculate the LO production rate both for a finite top quark mass, and in 
the limit that the top mass goes to infinity (the large irtt limit) . We will then argue 
that in the Higgs mass range under consideration, the finite mass calculation and 
the infinite-mass approximation give virtually the same result, whereas the QCD 
corrections give very large contributions to both the transverse momentum spectrum 
as well as the overall production rate. 

4.1 The Gluon Fusion Cross Section 

The gluon fusion process proceeds through a virtual top quark loop as represented 
by the Feynman diagram in figure f4.1| . As this process is well known in the literature 
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Figure 4.1: Feynman Diagram for Higgs Production through Gluon Fusion 



25| , we will just present the result 



The form factor A is 



where 



at 



a 



256nv' 



A(r) = r[l + (l-r)/(r)] 



(4.1) 



(4.2) 



/(r) = arcsin^ (y^) 
/(-) =-i(log(^ 



ITT 



for r < 1, 
for r > 1, 



(4.3) 
(4.4) 



and r = 4m^Qp/m|^. In the limit of a large quark mass, this form factor converges to 



a constant 



A{r ^oo) = --, (4.5) 
which results in the cross section in the limit of a large quark mass; 



cr r ^ oo 



57Qttv^ 



(4.6) 



A graph of the exact hadronic cross section is given in Figure [4.2| for a top 
mass of 174 GeV. In addition, the corresponding infinite mass limit is plotted on the 
same graph as a reference. Clearly, both graphs are qualitatively the same for an 
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Figure 4.2: Comparison of the infinite quark mass approximation with the mtop =174 
GeV calculation. 
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intermediate mass Higgs {niH ~ 100-200 GeV), which is the region that we will be 
considering. If we assume that the LHC will run at a luminosity of approximately 
~ 300 fb~^ yr^^ then we may compute the expected number of Higgs produced 
in one year of continuous running by multiplying this figure by the cross section. 
Looking at the data on this graph, we see that the LHC will produce approximately 
0.4 — 4 X 10^ Higgs bosons in a year of running, depending on the Higgs mass. We 
can conclude that Higgs production will be quite prolific at this machine. 

The close correspondence of the finite calculation with the infinite mass limit 
calculation seems to suggest that in the infinite quark mass limit, the dynamics of 
the loop freezes out, allowing us to replace the loop with an effective coupling. As 
long as the relevant energy scales in the problem do not become too large, we can be 
confident that the results of this calculation give a fairly good approximation to the 
true cross section. 

Even in the region where the energy scale is approximately the top quark mass, 
the infinite mass approximation is quite good. The Ofm^p] give additive corrections 
26| which behave as 



^ = oo) + ^ j (4.7) 

so that for a Higgs mass with a mass as large as 200 GeV, the exact leading-order 
expression differs from the large mtop approximation by ~ 15%. Since every exper- 
imental indicator tells us that the Higgs is very likely below this limit, we can with 
good faith apply this effective theory to the inclusive cross section. 

We may also note that the cross section calculated in the large mtop limit 
underestimates the Higgs mass over the entire range of plotted masses. As we argued 
earlier, it is highly improbable that the Higgs mass is outside of this range. We 
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therefore expect that the infinite mass effective theory will be quite good over the 



entire range of applicable masses. Moreover, recent work ^ seems to suggest 
that the ratio of the next-to-leading order calculation to the leading order calculation 
is more or less independent of the top quark mass. Therefore, one may compute the 
next-to- leading order corrections using the infinite top mass limit, divide that value 
by the leading order, infinite mass limit result, and apply the resulting enhancement 
to the leading order expression with confidence. 



4.2 Low Energy Effective Lagrangian 

Since the large mtop approximation is quite valid in our leading order calculation of 
the cross section, a systematic way to calculate the various Higgs observables using 
this approximation is desirable. We may therefore propose a low energy effective 
theory to do gluon fusion calculations []. 

This idea of freezing the top quark inside the loop, and replacing it with an 
effective coupling leads to the following effective Lagrangian, 

where the top-loop dynamics have been replaced by an effective coupling 0, In 
order to be valid, the relevant energy scales of the process must be less than the top 
mass, and even outside that neighborhood this approximation generally underesti- 
mates the cross section, so that one errs on the conservative side. We may derive the 

^ Aside from the original leading order calculation [^9| , The radiative corrections using this infinite 
mtop limit have also been computed p6| , Moreover, progress has recently been made on the 
NNLO calculation |p, |3|. 

^The effective coupling a^/STrw receives a finite renormalization at each order in perturbation 
theory. We represent this with a constant A, which is zero at leading order, there is a contribution 
from A at NLO which one must take care to include. 
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following Feynman rules from this effective low energy Lagrangian: 
The Higgs-gluon-gluon coupling is 



(4.9) 



the Higgs-three gluon coupling is 



H^7{pi,P2, Pz) = gAUc {{Pi - P2Y 9'" + (P2 - PzY r + {P3-Pir gn , (4.io) 



and the Higgs-four gluon coupling is 



(4.11) 



where A — (1 + ^) the effective coupling and fahc are the S'C/(3) group 

coefficients in the adjoint representation. These are represented graphically as: 




^abc 



a, jj I u, ty 



b, u 




a, n ^abcd a, a 



In the preceding expressions and diagrams, the gluon a has color a, momentum pi, 
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and helicity index yU, with analogous relations for the remaining particles. Note that 
these rules become those for QCD if one simply takes the limit pn —>■ and replaces 
the coupling A with 1. 



4.3 The Leading Order p_\_ spectrum in the large 

rritop hmit 

By far, the most significant contribution to Higgs production is through the giuonic 
coupling via a virtual top quark loop. Not only is the rate enhanced by the large 
top Yukawa coupling, it is also enhanced considerably by the large number of gluons 
(in comparison to quarks) inside the proton. Therefore, in computing the transverse 
momentum spectrum, we shall restrict our calculation to radiative corrections of the 
gg H cross section. In particular, those that produce a nontrivial transverse 
momentum will be considered. When we refer to the Leading Order or Tree Level 
contribution, we mean the lowest order in perturbation theory which produces this 
nontrivial contribution. 

At a proton-proton collider such as the LHC, by far the most dominant con- 
tribution is the process gg Hg. This is simply because the gluon distributions 
are much larger than the quark distributions at these energies. The quark anti-quark 
mechanism is also small, since the antiquark distributions are also very small. 

We would like to use our effective theory to calculate this process, but we 
have provided no arguments that the effective theory will be any good whatsoever 
in this markedly different calculation. Therefore, we will compare it to the finite top 
mass calculation done in reference [^|. Indeed, we will present the argument that 



the effective theory does quite nicely as long as the mass scales in the problem (in 
this case, mu and p±) are fairly low in comparison to the top quark mass. 
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4.3.1 Kinematics and Notation 



We begin by defining our kinematics, and establishing our notation for tlie calculation. 
As is the convention of parton model, we relate the hadron momenta. Pa, Pb to the 
parton momenta pb, Ph via momentum fractions, ^o, 

Pa = U Pa, Pb = ^b Pb- (4.12) 

We will sometimes use the convenient shorthand a — Pa, b — Pb ior the parton level 
momenta. It is important to note that the final state parton momentum will be called 
Q in this work. This differs from many similar calculations in which the momentum 
of the final-state massive particle is referred to as Q. We therefore have the following 
momentum conservation relation; 

Pa+Pb^PH + Q. (4.13) 

We shall define the parton Mandelstam invariants in the usual fashion: 

s = {Pa + Pbf, 

t^iPa- PhY, u^ipb- PhY- (4.14) 
If we define the hadronic Mandelstam variables analogously 

S^{Pa + Pb)\ T^{Pa-pH)\ U^iPb-pH?, (4.15) 
these are related to the partonic variables via 

t = mUl- Q + ibT, u = ml{l - ia) + ibU. (4.16) 
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Figure 4.3: Feynman diagrams for Higgs production in the large m± limit. 

If we choose to evaluate the Higgs momentum in the hadron-hadron center 
of momentum (COM) frame, as is usually the convention, we may re-express the 
hadronic variables as; 

^ — ^COMi 

T = m]j- m±EcoM exp {-yn) , U = m]j - m±EcoM exp (yn) ■ 

(4.17) 

Here, the rapidity of the Higgs in the COM frame is denoted yn, and ■m± denotes the 
transverse mass (= y'ml^ +Pi)- Therefore the partonic momenta can be expressed 
in terms of the traditional external variables and the momentum fractions; 



•5 — CaCb-^COM) 

t^m]j- ^am±EcoM exp (-yj?) , u^m]j- Cbm±EcoM exp (yn). 

(4.18) 



4.3.2 Analytic Procedure 

The Higgs p± spectrum has three contributing partonic processes at leading order, 
gg — > Hg, qq — > Hg, and gq Hq. These may be derived from the graphs in Figure 
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using the Feynman rules outlined above. We find the following expressions for 
the (helicity-averaged) matrix elements in = 4 — 2e dimensions 



^ ,9 24 4aM, ,s'^ + t'^ + + mrr , ^ , , 

2 16 ft^ + M^,, , „ 



AWI^ = 9(36) ^ — - - T • 



and 



M„a^„H^ = —— 1-e -2e — }. 4.21 

I 99-5^1 9(96)(l-e)7rt;2 \ -t ^ ^ -t j ^ ' 

Because of CP invariance, the process gq qH is identical to the expression for 
gq — * qH, and the matrix element for qq gH may be derived from the latter 
expression using crossing symmetry. The appropriate spin and color averaging factors 
must, of course, be substituted. These have been included in the preceding expression. 

The general expression for a 2 to 2 partonic differential cross section is given 

by 

, , , = MP. 4.22 

The hadronic cross section is related to the partonic cross section via the double 
convolution, as outlined in Chapter 3; 

d(^AB=J2^f d^afA^i{^a,fJ'F) [ d^bfB^ji^b, IJ'F)daij{^a, ^b) ■ (4.23) 

Using the expression for cross section found in the appendix D and the expression for 
the partonic cross section, Eq. ( ^^221) , we may derive the following expression for the 



differential cross section: 



dapp^Hj^ = EE fdUA^.i^l^F) rdUB^,{^b,l^F)5m T7^|M,,f . (4.24) 



dp± dy.. , J 
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Noting that the double convolution of the delta function may be written (also Ap- 
pendix D) 



where 



x± = — ^ — exp {±yH), (4.26) 



we arrive at our final expression for the differential cross section, 



dpi dyu Y i V-^-+ ^aS + T-Ml J,. ^,S + U-Ml^ 



xfA^,{Uf^F)fB-.A^b,f^F)j7—\Mi^\'. (4.27) 

lovrs 



Note that there is only one degree of freedom in the integration, as the momentum 
fractions are fixed by the on-shell condition in each term. 

4.3.3 Numerical Procedure and Results 

The convolution ( |4.27| ) must be evaluated numerically since the parton distribution 
functions are empirically determined. We shall use the CTEQ5L parton distributions 
| |34] |, and we shall choose our scale to be the transverse energy, m±. The number of 
(light) quarks is 5 throughout the calculation. 

The finite-mass calculation proceeds identically to the one outlined above. One 
simply replaces the effective theory expression for the matrix element with the exact 



result [|33|. Because the mass of the top quark is so large, we expect the large mtop 



limit to be quite good out to very large transverse momentum. 



Figure 4.4 displays a graph of the transverse momentum spectrum at zero 



rapidity for a Higgs of mass rrih = 115 GeV calculated using the full mtop dependence 
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Figure 4.4: Comparison of the exact leading-order spectrum with the effective 
theory calculation. 



(dotted line). Superimposed over this graph is the identical calculation performed 
with the simplification that the fermion mass is taken to be infinite. It is plain to see 
that the agreement is quite good out to the very large p± region, which is roughly 
about 300 GeV. There is really no reason to suspect that the behavior of any of the 
higher-order corrections will behave any differently than this. 

It is advantageous to break the infinite-mass limit calculation up into the 
various contributions from its partonic constituents, so that we may see which partonic 
graphs contribute the most to the cross section. A plot of the various contributions 
to the spectrum at leading order for a Higgs with a mass of 115 GeV is given in 



Figure 4.5. All of the terms exhibit the expected divergent behavior as p± goes to 



zero, with the exception of the s-channel piece, qq Hg. 

As can be plainly seen in the graph, the dominant contribution to the total 
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Figure 4.5: The Higgs p± spectrum at leading order 



cross section (solid black line) comes from the gluon initiated process, gg — > Hg 
(long dashes). This contribution is roughly an order of magnitude larger than the 
next largest contribution, qg ^ Hg (dashed line), in the small transverse momentum 
region. Even in the large p± region, the gluon initiated process is twice as large as the 
gluon-quark initiated one. The process qq — > Hg is completely negligible compared 
to these other two processes. 

To summarize, the gluon fusion process dominates at the LHC, and working 
with an effective field theory in which the heavy top loop is replaced with an effective 
vertex is a quite good approximation. Even for the inclusive cross section, the NLO 
pieces are found to produce enhancements which are much larger than the difference 
between the finite mtop calculation and this infinite limit, at least for Higgs masses 
in the expected mass range, 100-200 GeV. Since the differential cross sections at 
leading-order tend to be even more unreliable than the leading-order cross sections. 
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it is imperative that we include some of these next-to-leading order terms in order to 
improve our estimate for the differential cross section. 
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Chapter 5 

The Higgs Boson Transverse 
Momentum Spectrum: Gluonic 
Radiative Corrections 

In the previous chapter, we found that the dominant production mechanism for pro- 
ducing a Higgs boson at the LHC was that of gluon fusion, and that this process 
is well represented by an effective theory in which the massive quark in the loop is 
replaced by an effective coupling between the Higgs boson and two or more gluons. 

Moreover, the dominant partonic process is the gluon initiated gg — > Hg 
process, which is responsible for about 90% of the total cross section. This is primarily 
due to the large number of gluons in the protons compared to the quarks. 

As we have mentioned previously, leading-order calculations in QCD are at 
best an order of magnitude approximation to the actual process. This is reflected in 
a large scale dependence in the leading order cross section. If one wishes to have a 
more precise prediction for the cross section, one needs to carry out the calculation 
of the differential cross section to next-to-leading order in perturbative QCD. It is 
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hoped that the resuhing calculation will have less scale dependence, and therefore, 
less theoretical uncertainty. 

There are several different methods for doing this. One notable method is 



using a Monte-Carlo event generator to do the entire integration Since the 

matrix element contains soft and coUinear divergences, one must take care to extract 
the divergent pieces and integrate them in 4 — 2e dimensions in order to cancel the 



divergences in the virtual pieces. A calculation of this sort has been done ||3^ . 

Nonetheless, it is advantageous to do the entire calculation analytically for 
several reasons. First, one would like to calculate the process using conventional MS, 
and show that it agrees with the Monte Carlo calculation. Second, one needs an 
analytic formula in order to expand the matrix element around small transverse mo- 
mentum. This is necessary to calculate the B^"^^ coefficient in Collins-Soper-Sterman 
resummation scheme, and compare it to the result calculated in a different manner 



in reference To have a complete next to leading order calculation, good for 



all transverse momentum, the intermediate p_i_ calculation included herein should be 



combined with a resummed curve good in the small p± limit |^, ^ ^ ^ . 

In this chapter, we will be computing the dominant piece of the NLO correc- 
tions, the partonic process represented by gluonic radiation from the leading order 
process, gg —>■ Hg. This, we will argue, dominates the NLO cross section in the 
intermediate p± region with which we are concerned. 

To do this analytically, we will use the procedure outlined in Chapter 3. Since, 
in these corrections, we will be required to renormalize ultraviolet divergences and 
factorize collinear divergences, we will have to carefully compute the matrix elements 
in (i = 4 — 2e dimensions. The parton distributions that we will be using are calcu- 
lated from deep inelastic scattering data using the MS factorization/renormalization 
scheme, so we must use the same scheme in order to obtain a meaningful result. Once 
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we have a finite expression, the convolutions with the structure functions may be eval- 
uated numerically to obtain the cross section as a function of the various kinematic 
parameters. 



5.1 Notational Conventions 

The notations that we will be using are in complete analogy to those of the previous 
section. The 4-momenta of the hadrons shall be taken to be Pa and P5, and the 
4-momenta of the partons shall be denoted Pa, Pb (= ^a,b Pa,b) ■ The Higgs boson 
momentum shall be denoted H and the mass m^f, the transverse momentum p^, and 
the rapidity yn- Therefore, the hadronic (5", T, U) and partonic (s, u) Mandelstam 
variables have precisely the same form as the previous chapter. 

It is also convenient to define crossed momenta, pi 2 = —Pa,b, and to label the 
final state partons P3,P4- This allows us to express the invariants of the 2^3 process 
in terms of generalized Mandelstam variables. 



The invariant mass of the final state gluons and the Higgs invariants t — {a — HY, 
u — {b — Hy are related to these via 



Sij = {Pi+Pjf, 



Sijk = {Pi+Pj+Pkf- 



(5.1) 



t — S234 — 



(5.2) 



and 



u = S 



134 — 



(pi +P3+Pif- 



(5.3) 
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One may then show that the Higgs mass M is related to the remaining variables by 

^ s + t + u-Q^, (5.4) 



and that the transverse momentum p± of the Higgs may be expressed as 



= '"-Q'"'^'. (5.6) 



We will find it convenient to define the parameters Q±, v as follows: 



V = QVQi- (5.7) 

The physical interpretation of Qj, is that it is the transverse mass of the gluon pair. 
Finally, we will also use St and 5'^, which are defined as 

Su = -5i3 + -Si4 = - Q', (5.8) 

St = -S23 + ^24 = t-g', (5.9) 

and we will occasionally use the parameter S^^ 

^i/ = ^123 + ^124 = + g2 _ 5_ (5.10) 
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5.2 Cross Section at NLO 



To next-to-leading order, the two-to-two (2 — > 2) differential cross section is: 

dcTNLO = daf}^2 + o^4^2 + dafX^. (5.11) 

Here, d(T2^2 is the tree-level contribution, (i(T2^2 is the expression for the 1-loop 
virtual corrections, and daf}^^ is radiative corrections to the leading order process. 
If we include a coUinear counterterm d(T^2 ^^<i ^ UV counterterm d(T^2) then the 
cross section may be written as 



da 



dpi dyH 



= / j d^af{^a:l^F)d^bf{^a:l^F)s X 



' drr^"'^ dcT^^^ dcT^^ \ dcT^^^ dcj'^P 



Since these are to be computed using the conventional dimensional regular- 
ization, we must be careful to include all e dependence in the latter two terms. Note 
that the expression for the 2 — > 2 differential cross section becomes 



dt \ PI J 167rs2r(l-e) 



(5.13) 



The expressions from the other two 2 — > 2 terms (the virtual corrections and UV 
counterterm) are identical to this one with the exception that one must substitute 
the requisite matrix elements. 

The contribution coming from the radiative corrections (those involving 5 par- 
ticles) may be written 



dtdu \P\ ) ) r(l-2e) (27r)(47r)2 
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X 2 r de{sm{e)y-^' r d(l>{sm{cl)))-^'l-^. (5.14) 

Jo Jo iOTTS^ 

The contribution from the Altarelli-Parisi term is 

daj^s ^ fATTiiY as {AttY 1 
dtdu \pI) 2 7rer(l-e)2l67rs2 

/ -^P{zMza{s + t-M^)+u)\M{Uza,^b)\' + {t^u). (5.15) 

Jo Za 

In the preceding expressions, we must note that we have suppressed the sum over 
parton species for clarity. 

The angles ^, in the expression for the 2 — > 3 cross section will be the polar 
angles of the final state partons in the rest frame. The z axis in this frame will 
be defined by either initial state parton or the Higgs boson, whichever is convenient. 
See Appendix B for a further discussion on the explicit evaluation of these angular 
integrals. Once these integrations are calculated, we need only concern ourselves with 
the end result. 



5.3 Angular Integration of Real Emission Piece 

To calculate the real contribution, we must integrate the two angular degrees of 
freedom occuring in the matrix element. While there are a large number of terms to 
integrate, the calculation lends itself to an algorithm very nicely. The basic procedure 
for the angular integration is the following: 

1. Partial fraction the amplitudes using the relationship amongst the various Man- 
delstam invariants. This will reduce all of the integrals to the form (const) x 
^ijk^mn (const) X SfjS^^, for which standard formulae exist. Here, const 
is some collection of kinematic invariants that does not depend on the angular 
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variables 9 and 0. 



2. Extract the terms proportional to 1/(1 + e) and verify that they cancel amongst 
themselves. This verifies that there are no quadratically divergent pieces in the 
amplitude (which is true in general for QCD). 

3. Extract the universal — 1/e piece and verify that it cancels algebraically against 
the coUinear counterterm. 

4. Regularize the soft divergences (Q^ — > 0), and extract the universal soft pole. 
This will cancel against the soft pole in the virtual piece. 

Once this is done, the real matrix element is integrable over the convolution space 
outlined in Appendix D. The final integral must, of course, be done numerically. 

As discussed in chapter 3, the helicity amplitudes obtained in the literature are 
computed using gluons with two polarizations. To use MS, we must have expressions 
for the matrix elements with 2 — 2e gluon polarizations. Therefore, in order to use 
MS, we will need to calculate the 0[e\ corrections to these expressions due to higher 
dimensional gluon spin effects. This calculation is much simpler than the original 
matrix element calculation, since we need only to calculate the matrix element in the 
coUinear limit. We will explain this shortly. 

5.3.1 gg Hgg 

We now restrict ourselves to the process gg Hgg. To begin, we first break the 
cross section into two sections, a piece which represents the sum over the (+ + ++) 
and ( helicity configurations, and a piece which represents a sum over the 
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odd-helicity components, ( — Let 



dtdu \ dtdu I , V dtdu 



+ ^^^^f^ = (5-16) 



with 



dagg^Hgg ] ^ f 

^ dtdu ^ ' \47ri^y J IGttV 

2 (|M(l+2+3+4+)|2 + |M(l-2+3+4-)|V 

|M(l+2-3+4-)|2 + |M(l+2+3-4-)|2) , (5.17) 



'dagg^Hgg\ _ r(l-e)(^V f^X 



and 



y dtdu \47r//^y J IGtt^s^ 

2 (|M(l-2+3+4+)|^ + 3 cyclic spin perms) . (5.18) 

In these expressions, the remaining parts of the 2 — > 3 phase space d^ may be 
succinctly written 

87r2r(l-2e) 1, J ' ^^^ ^^ 

where dQ^^'^ is the 4 — 2e angular integration measure, 

dfl^'^ = — r deisin (6))^-^' r d(j)(sm U))-^', (5.20) 
2tt Jo Jo 



which is defined in Appendix B. 
First, we define 



A^^ f as \ ( A-n iJi\ ( A-K ji' 



2567rs2 Ar2 _ 1 \^37r't;2y \^ g2 y \p\ ) 1(1 -2e)' 



(5.21) 
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Then, the angular integral of the first contribution in (|5.16| ) yields 



'da, 



dtdu 



, ,1\ M^ + s^ + St + Sf + Q^ 



;i + z/)(i-z/)-^ + (i-z/)z/-^(i 



+ 



St + Sf / 67e 1-u 



6 

y + (l~iyY-l)\ 



+ 4 



wi n\ e 17 

+ ln l-z/ + — + — 

1 — z/ / 3z/ 3 



(5.22) 



Noteworthy are the terms that behave like '"["L^'^'' ^^(^ p"? which diverge as p± 0, 
as well as the universal ^ pole. 

For the (—,+,+,+) terms, we write the results as a sum of the universal 
divergent pieces and an analytic piece which is finite as Q and p± go to zero. We 
find from the symmetries of the matrix element that there are in actuality only 4 
independent terms, and that the sum over the odd-helicity components may written 



colors, (—,+,+,+) 



4 I (^m(l-, 2+, 3+, 4+) + ^m(l-, 4+, 2+, 3+) + {u ^ t)J 
+2 m(3-,4+,l+,2+) + m(3-,2+,4+,l+)}. (5.23) 



Using this relationship, we may write 



'da, 



j = 2 a, {[2 Q{1, 2, 3, 4) + 2 Q{3, 4, 1, 2) + 4, 2, 3) + fi(3, 2, 4, 1)] 



+ {u ^ t)} 



(5.24) 



where the f2's are defined as the angular integrations of the helicity subamplitudes, 
e.g. 

f](l,2,3,4) = J ^|m(l+,2-,3-,4-)p, (5.25) 
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with similar relationships for the other subamplitudes. 

The angular integration over the (1-, 2+, 3+ 4+) hehcity configuration may 
be written as a sum over two terms, a term that is finite as Q, e — > which we shall 
call r2o(l, 2, 3, 4), and a term which contains the pole structure, which we shall call 
Qdi^.(l>2,3,4); i.e.. 



2, 3, 4) = no{l, 2, 3, 4) + (1, 2, 3, 4). 



(5.26) 



The second term in this equation is 



Qdi^.(l,2,3,4) 



1 



1 



2MH^ 



es Qlu{M'^ -t){M^ -u) 



-l-e 



u{M^ - u){M^ - t) 
(5.27) 



and the term flo{l, 2, 3, 4) may be found in Appendix C. The matrix element with the 
same color ordering, but with the final-state gluons having the opposing helicities, 
integrates to 



Odi..(3,4,l,2) = 



su 



2Q^sM\M^ - uf 
(M2 - t)tu 



[ t ) ^[it-Q')t^ 
e (M2 - t)stuQl ^ ' 



11 67e\ M« + e 



(5.28) 



The term with the complicated pole structure (1324) yields 



Qd«..(l,3,2,4) 



1 



u 



-2- 



+ffdy+,4(i_2fi) 

\ u tu 



2 2(1 - i/)-i-W^i3 



u{M'^ -u){M'^-t)\ es Qlu{M^ - t) (M^ - u) 
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for the term where the initial state particles have the opposing helicities, and 



ndi..(3,2,4,l) 



su V ^ f M\Q'-t) y ^ ( sQVi. 



u-Q'^J \ t J \{t-Q^)t^ 



(M2 - t)tu tu 



e{M^ -t)stuQY ' e sQltu ^^ 6^ 

+ {u ^ t)} (5.30) 



for the terms that have the final-state gluons with opposing helicities. Again, the 
expressions for the Qo's may be found in Appendix C. 

In using the helicity amplitudes to calculate the amplitude squared, we have 
neglected terms of 0[e] which appear naturally from the sum over spins. If these 0[e] 
spin pieces are proportional to a coUinear singularity, they will contribute to the finite 
pieces left over from the coUinear cancellation, and hence, the cross section. Therefore, 
it is necessary to calculate these terms in order to complete the next-to-leading order 
calculation. 

While it may seem formidable to calculate these terms, one realizes that one 
only needs the terms proportional to coUinear poles in the matrix element, since the 
remaining terms vanish when the limit e — is taken. The correction term may be 
calculated in this manner. For the (1234) color ordering, we find the correction term 
to be 



^M^r.-MAr,'j^^[ —^^-' ) ^^^j .^|. (5.31) 

In this expression, everything save the (5'i45'23 — 313824)'^ / S^^ is evaluated in the limit 



66 



5'34 0. The remaining divergent pieces (as S23 and S'14 0) may be obtained 
from this one by cychcally permuting the momenta in the above expression for the 
amphtude, and taking the subsequent hmit. Upon integrating this term, we find the 
following correction to the (1234) piece, 



where 




(5.32) 

(5.33) 
(5.34) 



Ts^T2{u^t). (5.35) 

The integral over the (1342) ordering yields precisely the same result. The (1423) 
piece simply gives 

E / (^) l\M^.\l2s - (-26 (Afg'^^^) m + 2T,) . (5.36) 

Note that, when added to the previous expression, all of the extraneous pole behavior 
of the form e/Q^ is eliminated. 
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5.4 Extraction of Soft Divergences and Collinear 



Renormalization 

The expression for the real piece contains coUinear singularities, which must be ab- 
sorbed into the PDF's. The expression for the splitting function is 



at AP IbiTS Jo z elTxl[\ — Ze\ 



1 dz I a^Pa'^aiz) 



|Mi5r(^ea',6) {^n)-^Siz{u - Q') - u). (5.37) 



IGnsJo z e27rr[l-2e] 
Using the S function to evaluate the z integrals gives the following: 

da _ (Att)-' lasPb'^bizt) iMi'^hHUzt^b') 



s 



dtAP 16ns e27rr[l-2e] -t 

IQns e27rr[l-2el -u ' ^ ^ ^ 



where the parton momentum fractions are defined via the delta functions to be 



t,u 

zt,u = 5.39 

t,u — Q'^ 



For the process in question, the matrix elements above become 



,2 x4 



while the splitting kernel Pg_,g{zt) is 



P,^,{z,) = ^^^^]^ ''^'^-^'^^ + PoS{l - z,). (5.41) 

Zt[i~ — Zt), 



The second term in (|5.38|) is related to this one hj u ^ t. 
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The entire pole behavior of the matrix element may be expressed as 



da 

dt du 



(J, 



-t 



— u 



1 

e 
1 

~2 
+Reg 



1 



1 



11 
~6 



6 J J Zts'^p'j_Q^ 

Q7e\ + jzts)^ + + jztspl/t)' 
18 J zIs^pIQ^ 



(5.42) 



where "Reg." represents all of the remaining nonsingular (as —>■ 0) pieces of the 
cross section. The functions Pg^g are defined to be 



:i+z" + (i-^) 

zd-z) 



(5.43) 



This term, when multiplied by the Q"^*^ implicit in a^, contains both soft and collinear 
poles. We must take care to extract the poles from this term, and cancel them 
accordingly. In addition, the term in the second line is also poorly behaved as 0, 
although it is free of any collinear divergences. 



In order to extract the infrared pole from the expression ( 5.42| ), so that we may 
cancel the remaining (collinear) pole against the Altarelli-Parisi expression ( p.38| ), we 
make use of the identity 



l+ne 



1 ^ 

nt 



+ 



Zt 



l+ne 



Zt 



ne 



'ln(l 



Zt 



This allows us to write the pole term as 



da 

dt du 



a. 



TT 

12 



-us: 



+ 



A + 111' 

eN^ 2e 6 , 



MS + ^4 + ^4 + ^4 
' ^1 
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1 Pg^,{zt) + (zts)' + t' + {ztspl/tY 



Nee 



+ 



ZtS^Pl 



'log(l - Zt) 



log 



ztn^ 



^- Zt J ^ \ t J \1~ Zt 

1 P^%{zt) + {ztsY + t' + {ztspl/ty 



N,e -t 
' logjl -Zt) \ 
. ^-zt I 



ZtS'^p']_ 



- log 



zjQX 



l-zt 



X 



M'' + {ztsY + t'' + {ztspl/tf 

ZtS^Pl 

111 f 1 \ + (zts)^ + t^ + {ztspl/tY 



2 6 \l-zt 
+Reg., 



zts^pli-t) 



(5.44) 



where Pg^g{z) is the gluon splitting kernel, and 



l + z^ + (l-z) 



(5.45) 



The second line now cancels identically against the gluonic contribution to the Altarelli- 
Parisi counterterm; equation (|5.38|) . The remainder is finite with the exception of the 
infrared pole proportional to the delta function. It should be noted that while this 
expression is now well behaved for any invariant mass Q^, it is poorly behaved as 
P±^0. 



5.5 Virtual Contribution 



The second term that contributes to the NLO cross section is the one-loop virtual 



correction [^]. This expression contains both ultraviolet and soft divergences. The ul- 
traviolet divergence cancels against the charge renormalization counterterm, whereas 
the soft divergence cancels against the soft divergence in the real contribution. 

70 



The complete 1-loop correction, including the ultraviolet renormalization piece, 
may be written as follows: 



da 
dtdu 



. (, ir^ ,Af» + a" + + 



(5.46) 



In this expression, is the universal singular piece, which has been renormalized as 
described in Chapter 3 and Appendix B. Explicitly, it is 



+ 



-u 



71^ 11 3/?o 



6 eN, 



+ log 
-log 



m 



2 \ 2 
H 



m 



+ log 



log 



m 



H 



m'\j — t ^ 



log 



H 



my — u , 



-t 



H , 



m 



log 



+2Li2 1 - 



m 



H 



2Li2 



m 



log 



-u 



H , 



m 



log 



H , 



m 



log 



-t 



H , 



m 



H , 



m 



H 



mjj — t , 



+ 2Li2 



H 



my — u 



(5.47) 



In this expression, we have rewritten the dilogarithm functions Li2 using several 
transformation formulae. The above expression is completely real over the physical 
range of the parameters. 



5.6 Cancellation of Soft Singularities and 

Renormalization of Ultraviolet Singularities 



When the contribution from the virtual matrix element ( |5.47 ) is combined with the 



term proportional to 5((5^) from ( ^.42 ), we find that all of the singular terms cancel. 



up to terms of order n/, which we are neglecting. Therefore the finite remainder of 
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the cancellation of the soft singularities Ureg is 



71^ 67 11 1, 11 / 



3 18 Nc 2 °\tj 12 \spij 
+ log +log^^ +log' ^ 



s j \mjj — tj \mjj — u J 

-log (^) log ) - log ( ^) log ( ^ ) - log (^) log ) 

+0[nj], (5.48) 



where we have set nj equal to zero. To complete the cancellation of the soft poles, the 
soft singularities resulting from the process gg — »• Hqq must be also included. This 
piece contributes the remaining nj/6 piece required to cancel the (complete) Pq pole. 

For the present calculation, we will assume that these fermionic pieces do 
cancel, and that the finite remainder contributes very little to the overall cross section. 
With this simplification, we take f3o Ncll/Q, and the pole cancellation is complete. 
While this will probably give us a very good approximation to the cross section, the 
gluon initiated process gg Hqq is probably of the same order of magnitude as the 



pure gluonic process |H3|, and must eventually be included in a complete analysis. 



5.7 Numerical Procedure and Results 

In order to generate a transverse-momentum spectrum that may be observed in na- 
ture, we must numerically evaluate the convolutions of the preceding quantities with 
the parton distribution functions. We choose to use the CTEQ5 structure functions 



34 . The convolutions of the matrix elements with these structure functions must 



be evaluated numerically. We use a simple Gaussian integration routine, as found in 
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Numerical Recipes for C [Q. 

With this code, we may generate the differential cross section for any transverse 
momentum or rapidity (assuming of course that these are kinematically allowed). The 
process is to be plotted as observed at the LHC, a proton-proton collider running at 
a COM energy of = 14 TeV. Unless otherwise stated, we shall use a Higgs mass 
of 115 GeV in our evaluations. 

The leading order cross section, when shown, was computed using the CTEQ5L 
distributions, together with the leading order definition of Ug. The next-to- leading 
order cross section will always be computed using the NLO definitions of the CTEQ 
PDF's, namely the CTEQ5M1 distributions. In particular, we have taken care to 
define the K factor as it is traditionally defined, as the ratio of the NLO cross section 
(with NLO PDF's and as) to the LO cross section (with LO PDF's and as)- Also, 
all contributions involving quarks have been ignored in this calculation Q We do 
however, keep n/ = 5 in the parton distributions and the strong couphng as- 

We shall begin with a discussion of the transverse momentum spectrum, fo- 
cusing on how it differs from the leading order behavior. Then, the scale dependence 
is addressed. We close by looking at the dependence on rapidity and the behavior of 
the cross section at small transverse momentum. 



5.7.1 Spectrum and the K Factor 

Theoretically, the dramatic scale dependence of the cross section at LO should be 

reduced significantly at NLO. Figure [SH] shows the cross section at NLO superimposed 

over the cross section at leading order over a wide range of rapidity. We have chosen 

to plot bands in this case in which the scale has been varied from .5m_L — > 2m±. 

The leading order graph (light grey) clearly varies more over the chosen scale range 
^For the sake of consistency, we use only the gluon-initiated graphs in the LO piece. 
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Transverse Momentum, GeV 



Figure 5.1: Variation of the p± spectrum with scale 

than the next-to-leading order graph. Again, this improvement is exactly what one 
expects for an NLO calculation, as the scale dependence represents the theoretical 
uncertainty in the calculation. As we improve our calculation order-by-order in QCD, 
our dependence on this artificial parameter should decrease. 

The graph of Figure |5.2| displays the dependence of the next-to-leading-order 
cross section on scale (dashed line) compared with the scale dependence of the leading- 
order calculation (solid line). We display the ratio of the cross section at a scale 
fi = Am_|_, normalized to the cross section for A = 1. The quantity plotted, then, is 
basically the fractional uncertainty in each distribution which is purely due to theory. 

As can be plainly seen in the graph, the leading order cross section increases 
monotonically as A ^ 0, whereas the NLO cross section appears to be approaching 
a peak around da/dao = 1.4. Eventually, the NLO graph will turn over and begin 
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Figure 5.2: Scale dependence at p± = 50GeV 

to drop to zero as the scale approaches Aqqd- However, it is very clear the NLO 
calculation represents a dramatic improvement in the theoretical uncertainty due to 
scale dependence. For instance, for a typical range .5 < A < 2, 0, the leading order 
cross section varies from 75% to 140% of the A = 1 value, whereas the NLO cross 
section only varies from 85% to 115%. 



In looking at Figure |5.1| , it is clear that the enhancement over the leading order 
approximation is very significant. The graph of the NLO cross section divided by the 
LO cross section (called the K Factor) is, perhaps more interesting, since it entails 
explicitly the magnitude of these corrections relative to the LO approximation. This 



is displayed in Figure 5.3 



The K Factor as a function of p± is fairly flat down to about 25 GeV, where 
it begins a dramatic decline. This is to be expected, since for small p±, the NLO p± 
spectrum behaves roughly like log^ {p±/'mH)/p'±, which diverges to negative infinity, 
while the leading order expression goes to positive infinity for small p±. 
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We should note that our expression disagrees substantially with the previous 
Monte Carlo calculation Whereas their K factor is about 1.5, ours is roughly 
1.9. The most likely reason for this discrepancy is that we have not included the 
quark graphs in our calculation, whereas the Monte Carlo calculation has these pieces 
included. It should be noted that our K factor is computed a.t yn = 0, whereas the 
Monte Carlo calculation integrates over y^- As we shall soon see, the K factor appears 
to exhibit a minimum at yn = 0, and therefore the integral over rapidity is unlikely 
to be the source of the disagreement. Were we to compute the integrated K factor, 
it would most likely disagree by an even wider margin. However, before claiming a 
discrepancy with the calculation in reference [Q, we must include the quark pieces 
in our calculation, or at least establish that their contribution is non-negative. 
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Figure 5.4: NLO Higgs p± spectrum at the LHC (v^ = 14 TeV) for various values 
of the Higgs rapidity, yn 

5.7.2 Rapidity Spectrum 



Kinematically, it is simply less probable to produce a massive particle at large rapidity 
unless the transverse momentum is small. Therefore, the cross section should fall off 
rapidly for large rapidity. Figure |5.4| is a graph of the complete p± spectrum at next- 
to-leading order for various rapidities, i/h = 0, 1,2, which are the solid, dashed, and 
dotted lines respectively. The graphs indeed decrease rapidly as the absolute value of 
the rapidity is increased. 



Figure ^]5| is a plot of the complete next-to-leading order rapidity spectrum 
at the LHC for p± values of 50 GeV and 100 GeV. This graph displays the classic 
peak at i/h = 0, and rapidly drops to zero as the edges of accessible phase space is 
approached. That is, the maximum rapidity of a Higgs Boson is governed by 



(p_L + m±) exp IhhI < y/S 



(5.49) 
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Figure 5.5: Rapidity Spectrum at the LHC (v^ = 14 TeV). The dotted hues are the 
leading order calc, and the sohd hues are the NLO calc. 

beyond which it is kinematically forbidden to produce a Higgs. 

The rapidity dependence of the K factor is plotted in Figure |5.6| . The first 
thing to note is that the graph is very flat as a function of rapidity, indicating that the 
phase space limitations more or less govern the shape of this graph. Also, note the 
minimum that we mentioned earlier at yn = 0, indicating that an integral over the 
rapidity will further enhance our expression for the K factor. Were we to integrate 
over i/H, this would most likely increase our K factor slightly. Again, this is interesting 
due to our disagreement with p6| . 



5.7.3 Small pj_ Behavior 

The transverse momentum spectrum calculation at any fixed order in perturbation 
theory is poorly behaved for very small p±, as it is dominated by logarithms of the 
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Figure 5.6: Rapidity dependence of the K factor (v^ = 14 TeV) 

transverse momentum divided by the Higgs mass in this hmit. At each order in 
perturbation theory, the coefficient of the leading divergence changes sign, therefore 
rendering the predictive power of traditional perturbation theory useless in this limit. 

For our calculation, we expect to see exactly this. The cross section should 
display a peak in the small p± region, after which the cross section starts to diverge to 
negative infinity. We expect that our calculation will be a fairly good approximation 
beyond this peak. 

Our plot. Figure |5.7| , of the small p± region displays just this behavior (dotted 
line). Again we have superimposed the leading order calculation (solid line) as a 
reference. Our calculation displays a prominant peak around 2.1 GeV for our standard 
choice of the Higgs mass, 115 GeV. 

If one wishes to have a good prediction of the cross section in the small p± 
region, one may resum the large logarithmic behavior via the techniques developed 
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Figure 5.7: Small p± dependence of the cross section, yu = 0. 



by Collins, and Soper and Sterman The first order of business for a cal- 
culation of this nature would be to expand our analytic result, and extract the -B*-^-* 
resummation coefficient. This could then be used to verify the result of reference 
||37|| . In order to have a complete small transverse momentum calculation, one needs 
to analytically calculate the contributions from the process gg —>■ qqH in the small 
p± limit. 

This work is a major part of a larger calculation. Whereas we have the dom- 
inant part of the cross section, the gluonic piece, completed, we need the quark 
contributions to have a complete and reliable prediction. One also wishes to have 
a reliable calculation at small transverse momentum, which involves calculating the 
resummation coefficient mentioned above. Including the quarks is work in progress, 
as is the computation of the resummation coefficient from the asymptotic expansion 
of our analytic results. 
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Chapter 6 



Conclusions 



In this work we have seen how the Standard Model of particle physics, while close to 
being complete from a phenomenological standpoint, still lacks a crucial piece, the 
Higgs Boson. The advent of the LHC, which in effect will be a "Higgs factory" , will 
allow us to study this component of the SM in considerable detail. If we are to under- 
stand the mechanism behind electroweak symmetry breaking, we must understand 
the behavior of a pure Standard Model Higgs as much as possible. 

One of the most important properties to understand is the transverse momen- 
tum spectrum of the Higgs at the LHC. Primarily, one needs to know what to expect 
in a detector as far as the decay product energies, among other things. Since higher 
order corrections to p± spectra from QCD radiation tend to be quite significant, one 
may not consider the LO process to be reliable. 

As we have seen, the results herein have verified that perturbative corrections 
to the Higgs boson transverse momentum spectrum are very significant. The terms 
that we have calculated enhance the cross section by a factor of approximately 2 
in the large transverse momentum region. In addition, we find that the theoretical 
uncertainty due to varying the renormalization scale is far less than the leading order 
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result, indicating, that our calculation is more reliable. This result also agrees with 
the expected asymptotic behavior in the small p± region, and agrees qualitatively 
with the Monte Carlo calculation by DeFlorian et.al. ||3^. In order to precisely 
compare their the results with this calculation, the NLO quark terms, which have 
been neglected, must be included. 

In the next ten to twenty years, much will change as to our understanding of the 
fundamental building blocks of nature. In particular, many clues to the fundamental 
phenomenon behind electroweak symmetry breaking lie in wait to be uncovered. The 
discovery of the phenomenon behind the electroweak symmetry breaking mechanism 
will probably be the first significant increase in our understanding of particle physics 
since the development of the Standard Model in the 1970's. Discovering the Higgs, 
and understanding its properties, is the first crucial step. 
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Appendix A 
Matrix Elements 



In this appendix, we include the expressions for the (unintegrated) matrix elements 
for the 2 — > 2 and 2^3 gluonic parton processes, computed in the large quark 
mass limit. Moreover, we include the expressions for the 2 — > 2 process H —>■ gqq, so 
that the complete leading-order transverse momentum spectrum may be computed 
in the large mtop limit. H ggg and H gqq amplitudes have been known in the 
literature for a while p3|, as have the expressions entailing the process H gggg 



47 1 . Moreover, the amplitudes for the virtual corrections to the 2^2 partonic 



processes H ggg and H gqq may be found in reference . For simplicity, we 
give the matrix elements for all particles, except the Higgs, in the final state. The 
matrix elements relevant to our calculation can be obtained from these by crossing 
symmetry. 

Since we will be using color-ordered helicity amplitudes in our calculation 
we will begin by reviewing the standard notation and identities associated with the 
decomposition of QCD amphtudes into color-ordered, gauge-invariant subamplitudes. 
We will then present the matrix elements relevant to this calculation in their entirety. 
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A.l Helicity Amplitudes 



We begin by presenting the standard notation used in calculations of amplitudes, and 
list a series of identities useful in simplifying these same amplitudes ^. We introduce 
the following notation for (massless) Dirac spinors, 

u±{p) = \p±). (A.l) 

In addition, we shall use the following notation for spinor products, 

u-{p)u+{q) = {p - \q+) = {pq) , (A.2) 
u+{p)u^{q) = {p + \q-) = [pq]. (A.3) 

The momenta p, q are assumed to be massless, and therefore these may also be iden- 
tified with relativistically invariant chirality states. 

The polarization vectors of the gluons are parameterized in a dual basis ^ 



where p is the momentum of the gluon, and k is some arbitrary massless reference 
vector. 

Products of these spinors are found to satisfy the following identities (for 
massless p, q)\ 

{pq)[qp] = 2p ■ q (A.5) 
(pq)* = sgn{p ■ q)[qp] (A.6) 



^ Those who are new to hehcity spinor techniques will find references |l9| useful. 
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^^4^;^=b±)(p±l- (A.7) 



Prom this last relation, we see that we may derive the following identity, 

-P2)(-P2| +P3)(+P3| -Pa) ■ ■ ■ {-Pn\ = ^Tr{(l + 75)|^i|^2|^3|^4- • ^n] ■ 

(A.8) 

This last relation is particularly useful in computing absolute squares of matrix ele- 
ments, and converting them into traces over products of gamma matrices. 



A. 2 Real Contributions 

In order to calculate the spectrum to next-to-leading order, we need to write down 
all the separate partonic contributions to this process. They can be broken up into 
three separate terms; the leading order or tree- level process which consists of Higgs to 
3-partons , the virtual or 1-loop correction to this process, and the 4-parton process, 
which we refer to as the real corrections. The (unintegrated) expressions for these 
have been computed in works outlined below. 

One may compute the spin- and color-averaged matrix element squared by 
simply summing over colors, adding up the 2^ spin terms, then dividing by the 
number of spins and colors of the particles in the initial state. That is, if Na,b are the 
number of spins and Ca,h are the number of colors of the initial-state partons a and 
6, then 

|M(a,6;l,2,...,n)|^= \ ^ |M(a^^ 6^- . . . , n^")|^ (A.9) 

where the represent the helicity of particle i, and the sum is implicitly over the 
spins and colors, unless noted otherwise. Generally, the amplitudes for the process 
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involving initial state partons may be computed using crossing symmetry, and being 
careful to note that certain terms can pick up a sign via equation ( |A.6|) . 



A. 2.1 H ggg 

We shall begin by writing down our expressions for the Higgs production matrix 
elements. The matrix element with the full e-dependence was first computed by S. 
Dawson It is fairly straightforward to compute from the conventional Feynman 
rules (together with the non-conventional ones from our low energy Higgs model). 

There are 2 independent subprocesses that contribute to the leading-order 
expression. The H —>■ ggg piece is given by 

M(l+,2+,3+) = sAfabcM^ , . 

M(l-,2+,3+) = 9Afabc[2S]' ^ 

^ ^ V2[12][13] ^ ^ 

The remaining contributions may be obtained through parity and charge conjugation 
operations. If we sum this over spins and colors, we find 

E mH - 999)\' = ''^l^fl ~ {St, + St, + 5,^3 + M^) . (A.12) 

The remaining subprocess that we will consider for the leading-order transverse 
momentum spectrum is if — > qqg. The matrix element for this subprocess is given 
by 

M(1^2-,3-) = -^li|. (A.13) 
This may be also be easily summed over spins and colors 
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It should be noted that each of these amphtudes were computed using the conventional 
normalization for the SU(3) generators, Tr[TaTb] = ^Sab- 



A. 2. 2 H gggg 

In order to facilitate the calculation, it is beneficial to break up the subamplitude 
into color-ordered helicity subamplitudes 

M(1^^2^^3^^4^4) = ^ Tr[T"T^T"T'^]m(l^S2^2^3^«,4^^). (A.15) 

colors 

These amplitudes are cyclically symmetric in the colors, they are invariant under 
reversal of the color indices, and they satisfy the dual Ward Identities [|T9|. There are 
three independent helicity subamplitudes for this process: 



""^^■^^■^^■^^)= (12)(23)(34)(41) - f^-^^) 



124 '-'12 '-'14 'J 123 'J 12 '-'23 



(1-||^h|2-)^[34]^ [24] f {1-1^2-) 

Su^SuSu [12] (23) (34) [41] I (41) 



+534^^^§|^-[24]5234}, (A.17) 



(A.18) 



and 



mil- 2- 3+ 4+) = (A 19) 

^ ' ' ' ^ (12) (23) (34) (41) [12] [23] [34] [41]- ^ ^ ^ 

The remaining terms can be related to these by using the charge conjugation and 
time reversal transformations, and the dual Ward identities. 
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A.3 Virtual 



The expression for the virtual H ggg amphtude was calculated in reference |[42|| . 
We define the amplitudes in precisely the same way that we defined the leading-order 
pieces, in terms of helicity amplitudes. The resulting expressions are then related to 
the leading-order amplitudes via a coefficient. These are outlined in the following 
subsections. To construct the contribution to the amplitude squared, use the formula 



|jVf(o) + M«|2 = J2 + 2 E Re[MW (M«)*] + . . . 



(A.20) 



spins 



spins 



In this expression, M(°) and M(i) represent generic LO and 1-loop amplitudes. For 
clarity, the spin dependence has been suppressed in this equation. 



A.3.1 H ggg 



The expressions for the virtual H — > ggg subamplitudes are 



mW(1+,2+,3" 



NrU 



, N 'S'3iS'23 + S'3iS'i2 + 5'i25'23 



M«(l-,2+,3+) 



+ -{N,-nf)—^ 

<J '-'23 



(A.21) 



where the prefactor is 



rr 



r(l + er(l-e) 
rfl-2e) 



r(i-e) 

r(l-2e) 



(A.22) 
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and 



U 



In 



—Sl2 

-Sl2 



'23 



-Ml 



In 



-S- 



23 



-Ml 



- In 



-Ml 
—S31 
-S12 

l4l 



+ 



In 



2 

'S'si 



-Ml 



2Lio 1 



S 



12 



2Lio 1 



s. 



23 



Ml 



2Lio 1 



- In 

'Ml 



'23 



-Ml 



In 



-S: 



31 



-Ml 



(A.23) 



It should be noted that all of the singular behavior (including as p± —>■ 0) is contained 
within terms proportional to U. 

Since we are working with an effective theory, we obtain a finite renormaliza- 
tion at each order of perturbation theory, 



^ m(i) -f (A -f 3Sg) M(°) , 



(A.24) 



where A is the finite renormalization of the effective Hgg operator, at this order 
in perturbation theory, and 5g is the gauge-coupling counterterm. Using the MS 
subtraction scheme, the counterterm is 



1 

e 47r 



r(l + e) {AttY 



UN, 



1l 
3 



(A.25) 



and the finite renormalization coefficient is 



A = 11 



2 7r 



(A.26) 
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Appendix B 



Formulas for Angular Integration 



In this appendix, we give detailed formulas for the angular integration of all the terms 
appearing in the partial-fractioned matrix element. There are two different cases: 

• when the two invariants involve massless particles, and 

• when one of the invariants contains a mass term. 

Momentum conservation may be used to reduce all of the matrix elements into sums 
of terms of this form, multiplied by a coefficient which has no angular dependence. 



B.l General Expression for terms of the form 

'^13'^23 



We begin by evaluating integrals involving two massless invariants. There are two 
forms that occur, and these are related by a — > {1 — v) symmetry: 




(B.2) 
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We recall the explicit defintion of 

= (B.3) 

Those that involve the other combinations of the invariants are related to these via 
the freedom to relabel the final state parton indices, that is exchange 3 -^-^ 4. 

First, we evaluate the phase space in the rest frame. Then without los- 
ing any generality, we perform this integration assuming that the invariant 5'i3 is 
independent of the azimuthal angle The following definitions are useful: 

c^i3 = ^is/l-^ia + ^u), (B.4) 
= SxzIS^, (B.5) 

C^14 = -5i4/(5l3 + -Si4), (B.6) 

= Sx^lS^, (B.7) 

<^23 — 'S'23/('S'23 + 'S'24), (B-8) 

= S^^ISu (B.9) 

<^24 = 'S'24/ (5'23 + 5*24), (B.IO) 

= S^^ISt. (B.ll) 



Taking the gluon a to lie along the z axis, we find 



2 

and 



C^13 = J(l-cos(^i3)), (B.12) 



u;23 = ^(1 - cos(^o) cos(^) - sin(eo) sin(^) cos(0)). (B.13) 



91 



The integrals ( p.2|) become 



~m -n 
UJ12, UJ24: 



(B.14) 
(B.15) 



These may be evaluated in general in terms of hypergeometric functions The 
explicit expression (for arbitrary e) is 



r2(l - e) r(l - n - e)r(l - m - e) 
r(l - 2e) r(2-n-m-2e) 



Fi,2(m,n,l-e,l-z/). (B.16) 



One may easily prove that the relations for 3^4 yield precisely the same result. 



B.2 Expressions for terms of the form 



on cm 
'^13'^123 



For these expressions, a considerable simplification is attained if one rotates the pre- 
ceding frame so that the Higgs momentum lies in the direction of the z axis. In the 
rest frame, this is equivalent to a frame where a + b = z|a + b|. The rotation 
angle is then given by 

cos(x) = , ~^ , (B.17) 
and therefore the invariant 5123(4) is given by 



5'i23(4) = At VW+C^cosOh, (B.18) 



where we have defined 



A = ^{2s + S^ + St), (B.19) 
B = ^{Su + St cos Oo), (B.20) 
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C = ]^{Su + St COS 6q). 



(B.21) 



The kinematic invariants (S'la, etc. take on the requisite form 



5'i3 = — (1 — cos X cos 9h — sin x sin 9h cos 0) , 
5'23 = — (1 — COS X COS Oh — sin x sin Oh cos 0) . 



(B.22) 
(B.23) 
(B.24) 



where the other terms are related by S — > — S. 



We therefore seek to evaluate expressions of the form 



n,m / 



(B.25) 



All of the other angular form factors may be obtained from these by either exchanging 
B <-> —B ox u *r^t. Their explicit expressions are (to 0[e°]) 



"if 



1,2 



1] 



2,2 
H 



1 



1 



1 



log 



A + S'cosx"^ ' ' V2a7 A + 5'cosx 

\ Q0,1(,) + \ 

(A + 5'cosx)2 ' V2ay (A + 5'cosx)^ 

■ + cos x) ^ 2S'2 + 1AB' cos X 
log ^ + 



(A + 5' cos x)^ 

A2 - 



X 



(B.26) 
(B.27) 



- S'2 



A + cosx 

1 



- S'2 

/ A^^cosx + ^^' \ 
V(^ + S'cosx)^ y a 



2a2 



AS'cosx + S'^, (A + S'cosx)^ 
log ' 



(A + S'cosx)^ 



^2_S 



/2 



2(^^sinx)^ 
(A + S'cosx)^ 



(B.28) 

;b.29) 



(A + S'cosx)^ 



^ /3(S'cosx)^ + 25'cosx(^ + S'cosx)^l ^ 



0,1/ 



log 



(--l + f)"c.>S v)2 



2a 



(^ + S'cosx)^ 

l(^^smx)^ 
a2(>l + S' cosx)' 



B 



a2(A2-S'2)(A + S'cosx)^ 



(B.30) 
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1-1 
H 



1-2 
H 



2-2 



1-3 



I , A + B' 



2B 



1 



^2 _ 5/2 

a 



a 

B' 
a 

B' 



{B' + A cos xWh - cos X 



(5' + A cos x)^'if - cos 



^0,-1 



{B' + A cos x)^H - cos xf^^f" 



(5'2 _ ^2)^1,0 ^ ^ 



U,-l 



(S' + A cos x)^^if - cos x^H 



{asm xY 



1,0 



(S' + A cos x)^^^~^ - cos x^H 



-.0,-2 



+ 



3 (a sin x) 



2 r 



5/2 



- + An"£' + B'q) cosx 



5' 



{B' + A cos x)^^i} - cos x^H 



^1,-2 



+ 



3 (a sin x) 
2 



2 r 



- + AQ^-^ + S'(^) cosxn^° 



(B.31) 
(B.32) 
(B.33) 
(B.34) 
(B.35) 

(B.36) 
(B.37) 

(B.38) 
(B.39) 

(B.40) 
(B.41) 

, (B.42) 



where B' = :f^/WTC^. 
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Appendix C 



Results of Angular Integration 



C.l Expression for Regular Pieces of H ^ 9999 
Matrix element 



In this section, wc give tiic detailed results for the non-singular parts of complicated 

(H ) helicity configurations. Recall that the sum over spins of the H —>■ gggg 

amplitude may be written 

|m(l^i , , 4^3 ^ ^ p + |m(l^i , 4^^ , 2^^ , 3^* ) p } . (C.l) 

We write the expression as 
where 



m(l+,2+,3+,4+)|2 + |m(l+,3+,4+,2 



+ M2 



m(l+,4+,2+,3^ 



(C.3) 
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where 



m++__=|m(l+, 2+, 3-, 4-)|' + |m(l+, 3+, 4", 2-)|2 + |m(l+, 4+, 2", S'jl^ + 
|m(l+, 2-, 3+, 4-)|2 + |m(l+, 3", 4+, 2-)|' + |m(l+, 4", 2+, 3-)|' + 
|m(l+, 2", 3-, 4+)|2 + |m(l+, 3", 4", 2+)|2 + |m(l+, 4", 2", 3+)|l (C.4) 

The angular integrals of the preceding expressions are presented in their entirety in 
chapter 4. The remaining pieces, which we shall refer to as the odd helicity pieces, 
are 

m+++_=|m(l+, 2+, 3+, 4-)|' + |m(l+, 3+, 4+, 2-)|' + |m(l+, 4+, 2+, 3-)|' + 
|m(l+, 2+, 3-, 4+)|2 + |m(l+, 3+, 4", 2+) + |m(l+, 4+, 2", 3+) ^ + 
|m(l+, 2-, 3+, 4+)|' + |m(l+, 3", 4+, 2+) |' + |m(l+, 4', 2+, 3+) |' + 
|m(l+, 2", 3-, 4-)|2 + |m(l+, 3", 4", 2-)|2 + |m(l+, 4", 2", 3-)|2. (C.5) 

We define the normalized angular integrals of these expressions by 

n+++_=2 J dn^'^m+++_. (C.6) 

dfi^ is the angular measure in 4 — 2e dimensions from appendix B, 

/ dQW = ^ r de {sm{e)f-^'^ rd(/.(sin(0))(-2^). (C.7) 

J ZTT Jo Jo 

Using the various symmetries, as well as the fact that we may freely exchange 
the partons 3 and 4 (since this does not affect the outcome of the angular integration), 
we may rewrite this as 

n+++- = 2 2, 3, 4) + 2 Q(3, 4, 1, 2) + 3, 2, 4) + Q(3, 2, 4, 1) + (t ^ u), (C.8) 
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where 

Q{i,j,k,l)= J dQ%\m{i-,j+,k+,l+)\^. (C.9) 

Finally, we break these into a fairly simple singular term fl^{i,j, k, I) (given in 
chapter 4), and a nonsingular piece Q,o{i,j, k, I), which contains no e or singular- 
ities. We have chosen to include a number of (universal) nonsingular factors in the 
expressions for fl^{i,j, k, I). 



C.1.1 m(l-,2+,3+,4+) 
First, let us define 



'l--^. (C.IO) 



Then we may write the integral over the angular phase space as 



^0(1,2,3,4) = F^^\s,t,u) + F^'\s,t,u), (C.ll) 



where 



r \S,t,U) — 2 2 



ur su^ su su^ 



+ 



su {s + Su) Sh 

tMog( (^^^+g;j^;^^'^' ) tMog(^) (Q^u + sSh + uSh) 



+ . . ^ s ^ + 



log( 



sis + St)SH Q^su{s + Su)Sh 



s w 



2 



M2 Su^ ' \ Q^s^u^p^^ su PA_^ {s + Su) {s + St) 

M^ + {s + Su)^ 



su{s + Su)Sh s{s + St)SH 

{Q^u + sSh + uSh) \ 
Q'^ su (s + Su) Sh J ' 



(C.12) 
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and 



2M^ + (r(_i,i)(^,t)-r(_i,i)(t,n)) 
, 2M^(r(_,,i)(n,t) + r(_2,i)(t,n)) ^ 6M^r(_i,,)(^,t) 



The angular form factors in F^^\s,t,u) are 



r(-2,2)(^^, t) 



loglig) ^ 48 {Sl + spj-sQ'') {2M^Q^-u{t + u)) 



+ 



2^2.2 \ (52 + ^p2_^,g2)2_4^2g2.2 



-4sM2 



+ 



+ 



(C.14) 



^{-l,2){u,t) 



(C.15) 



r(-2,i)(^^, t) 



s2 (2 M^Q^-u{t + n))' 2 p| ^2 



V 



52 5j^-2sn5„ 2s2Q2p2_(5^^^_2s^)2 



(C.16) 



and 



^ " 5^ (-4M2 . + Sh') V -^M^s + Sn'- 



(C.17) 
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Finally, the coefficient is in this case 



-2t^ {t + uf 



C = 

" Q'^suSh^ Q^suSh^ 



(-4Q^ + 6Q2u) +2M^Q2 (4^t'^ -2tu-u^ + (* + «))] -(0.18) 



Note the existence of a "naive" Q"^ pole Qo(l, 2, 3, 4). It is superficial, since the limit 
as — >^ yields a finite result. 



C.1.2 m(3-,4+,l+,2+) 



Again, we write the integral over the angular phase space as a sum of logs, 
Qo(3,4,l,2) = Co + Cilogj^ + C,^^ 

The coefficients Cj in this expression are 

63g2 + i70s 8Q^s^ leg^s 27g^- 12g2s + i3s2 

" 3^ + sT^ + 

3g^-i8g^g + 52g=^ 4s^ _ 11 (g6-3g^g + 3g^g^-g=^) 

3g2st ^ 3^2 + sQ2u 3QHu ^ 

lig^S^ _ 38g2s2 _ 2g2s3 -9^2^^7452 _QQ2g2^^^3 

3St^ 3St^ ~3uSt^^ 3S^ ^ 3uSt^ ^ 
22 g6 + 66 gS + 129 g2 - 52 s3 -18 g^ s + 27 g2 s2 - 11 s3 

3g2s5t ^ 3g2M5t ^ 

23^ 2sSt 2sSt Sl_ 2St^ 11 Q^s^ UQ^_2Q^ 

3u'' ^ Q'^u^ 3^^ 3Q''u^ 3Sj ^ 3Sj ~3tSj^ 
2Q'^sSt s{3Q^ + 26s) 3Q^s'^-7s^ 3sSt St^ 

C~3 ' 2 ' Q J- c 2 ' Ei~2 ^ 9 c 2 "I" 

22g6 + 66g^S + 93g2s2_4s3 ^ _6Q4^^i5g2 <52_;^2g3 ^ 

3 g2 s 5*^ 3 g2 1 (S^i 

11 (g^ + 4gs + 6gS2^4g2g3^2g^) ^ (iig^ + 26g2g-2s2) ^ 
6 g2 s 5"* g2 s 5'„ 
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2{llQ^-s)St^ ^ 25Su ^ 24Q^sSu WQ^sS^ ^ 31sS^ ^ 24Q^Su^ ^ 
3 s s Sf^ Sf 

2St^ st^ st^ st^ 

nQ^S^ + 34Q^sS^-26s^S^ 2 {6Q^ S^- s S^- 13s^ S^) 
Q^sSt Q^st 

27g^^^-66Q^5^. + 26 5^,g^ 2 (ll - 135 

3st2 + 2>Q'^sSt 

2 (27 S^^-%s Su') 9Q^Sj + 26sS^' 16 {Q^ S^^ - s S^') 



+ 



and 



35^2 SQ^st st^ 

(C.20) 



Q'^st^u{M-^ -ty 

where the numerator A^i is 

TV, = (^[M'-Q'^)'' + s')s,'{s + S^)-Q''^{Qs'-S,{-2s' + 2sSt + S^') 

+Q^St^ {6s' + 2s^ (7St + SS^) 

+s' {34St^ + 34StS^-6Sj) 
+s {St + S^f (ll St" + 22 St Su- 5 Sj") 
+10 {St + S^) (3 St^ + 4StS^- Sj) 
+{St + Suf {St'' + 5StSu-Sj)) 
(2s^ + 6 St^ + St^ (31 s + 13 Su) 

+St' (aO + 20 sS^- 2 Sj^) 
-Su {2s + Su) {4s^ + Us^Su + llsSu^ + 3Sj) 
-2St {s + Su) (7s^ + 27s^ Su + 23sSj^ + 6Sj) 
+St'^ (l6s^-2Su (ll + 25 s Su + 9 Sj))) 
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St (4 St'' + St^ (29 s + 17 5„) + 2 St^ (2 s + 5„) (l2 + 9 s 5„ - Su^) 
+2 5^' (29s2 + 32s5„ + 9 5j) 
+2 5t {s + 5„) (3 - 5 5„ {s + 5„) (2 s + Su)) + 
(s + Su) (6 - S„ (12 + 28 Su + lQs Sj^ + 3 
-2g« (_25i^-5i=^ (7s + S„) + 5i2 (-3s2 + 8sS„ + 55„2^ + 
5i (5s + 3S„) (s2 + 2S„ (2s + ^„)) 
+ + (6s3 + 5„ {lQs^ + 3Su (As + Su)))) . 

(C.22) 

Finally, The coefficient in equation (|C.19|) is 

^ 10 ^ {{AP-Qr + ^*) 

Q'^stu 

As (g^-3Q^5 + 3g^) 2g, (-4 (g2_33) 3_6gg^ + g^2^) 

r-i ■ (^^.23) 

C.1.3 m(l-,3+,2+,4+) 



Recalling Sh = M"^ + s — Q^, we write the result of the integration over the angular 



phase space as 



_ 16M^gW^4MMog(i±2 

"oU,d,Z,4J — h 



^(M2-t) (g2-M) SlM^s 



(C.24) 
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The coefficients are 



-2 (m^ -s) t^u^ + 2M'^t^u'^ {2Q'^ + t + 2u) 



Sh s (M2 - u) u"^ Su 

+4M^° (2Q^ -2Q^u + 2Q'^u^ -u^) 

-4M^ [Q^ {t-5u) u + Q^ {t + 3u) 

-u^ {t + 3u) + Q'^u'^ (t + du)) 

{t + 2u) (2t^ + u^^ -u^ (2t^ + tu + 
+4M^u (q^ {t + u) + Q* if + tu-Au^) 

+2Q^u {2t^ + tu + ?>vf) -u^ {2t^ + 2tu + ?>u^^y^ (C.25) 



and 



Co = 



(M^ (2g^-2t-M) + t (-Q^ + t + M)) 
s (M2 - t) (g2 -t) {2M^ -t-u) u 



(C.26) 



C.1.4 m(3-,2+,4+, 1+) 



Tfie fast independent lieiicity ampiitude is m(3 , 2+, 4+, 1+). We write tlie finite 
remainder of tlie tfie azimutfial integration as 

Qo(3,2,4,l) = -Co + Cilog. 



2 " ' ^M252 



u 



wliere 

/ {M^-tf \ {2t + u) {t^ + tu + u^) 2 M"^ {QH + {7 1 - u) u^) 



^St (-M2 + m) y stM StU^ 
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stu^ stu^ 



SU Su 



(C.28) 



and finally, the coefficient Cq, which is symmetric under the interchange u t: 



Co = 5 + 



st^ su^ su^ 

(M^{t-uf + (t + u) {t + 3u) + MH {-t'^ + tu + Au'^)) 



st"^ 



+ ^ - + 



2 {3M^Q^ -2M^QH) M^Q^u^2M^Q^u 



stu t^St'^ tSt'^ 

-10 Q^u-SM^^u + Q^u'^ + AQ^tu^ + eu'^ 

6M^t{2M^ + t) 2t{Q^ + 2t)u 10 M"^ + AM^t + 



o 3 o 3 ' o ^ 

iJt i~^t 

2M^t{M^-u) 4:M^{M'^ + t) 2M^t M^{-8M^ + t)St 
^ suSt" tSt suSt ^ ? 

6Q^{M^-tf 6Q^ {M^ -t) {Q^ + t) 4M^-8M^Q^ + g^ 

"I ^~3 I 

iJu iJu >Ju 

^2 g2 2 (2 M2 - 5 t ^ t2 

t Su Su Su 

2 {MH {t-3u)+t {Q^ + 2t) u + M^ (-2t + u)) 

+ I^u ■ ^ ' 



103 



Appendix D 



Formulas for Cross Sections 



Here we provide standard formulas for the evaluation of 2 ^ 2 and 2^3 processes 
for 2 partons in the initial state, 1 massive particle in the final state, and 1 or 2 
massless particles in the final state. We also hst convenient variable transformations 
for the convolutions, which will allow us to directly evaluate the "+" functions. 



D.l Expression for a 2 ^ 2 Cross Section 



The 2 — > 2 contribution to the p± may be written 



da 



LO 



1 /■! /■! drr 



dp± dy_ 
where 

g2 ^ s + i + ii - (D.2) 

and 

dt 47rs2' ' ' ^ ' 
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where we have suppressed the parton indices for clarity. 

The expression for the NLO cross section may be defined hkewise: 

dp± dyu S 

where 



JO Jo at au 



dan 1 /■ dVt I 2 



dt du Ans^ J 27t 



\M\'. (D.5) 



D.2 Evaluating the Convolutions 

First, we must evaluate the expression 

/' dU t d^h e{Q^) = r d^a t 0{Ub - XT^a - xu^b + t), (D.6) 
Jo Jo Jo Jo 



where, as before, Xt = Xu = — and r = Evaluating these in the 

hadron-hadron COM frame, we find that these expressions become 



Substituting expressions ( p.7|) for Q^, we find 



which has the implication that 



(eae-^--^)fee---^)>§. (D.9) 
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We may break this up into four different phase space regions, each defined by {^a,b^^^" — 
m^l \fS) greater than or less than \fS. One case of the four is automatically for- 
bidden because it implies p_L < 0. After evaluating the 6'-function for each term, we 
may recombine two of the regions into one region, leaving us with 



/' C di, e{Q^) = [' d^a t + /' di, r d^a. (D.IO) 



The lower limits of the internal integrals in this expression are solutions to the equa- 
tion = for each of the variables. Explicitly, these are 

C ic N (a^T - T)ia - r{l - ia) 

t,b{U> = — 



XT - 6 



whereas the outer limits x± are 



x^ = exp{±yH)r^^^]. (D.12) 



It is quite trivial to convert these terms back into the variable. If this is done, the 



expression for the phase space in reference BO, EOl] is recovered identically: 



Jo Jx+ q,a — Xu Jo \ J J^- — Xt Jo \ J 



In this formula, the limits of integration Ai are 

Ai = W^-xt)-xu + t, (D.14) 
A = x+ - xt) - xuib + r. (D.15) 

Note that it is quite trivial to evaluate the 5{Q'^) of the 2^2 phase space using this 
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form of the convolution. 



D.2.1 Evaluating Terms With "+" Functions 

In the regularization of the soft singularities, terms of the form f{zt^u)+ arise, and it 
is beneficial to change variables from d{Q'^) to dzt^u in order to easily evaluate these 
"+" functions. Recall the definition of Zt^u'- 

Now, consider Zf (a completely analogous analysis will follow for z^). It is beneficial to 
first convert to the phase space (|D.ltj|) as an intermediate step. We wish to evaluate 



the Jacobian for the change of variables. This is summarily expressed as 



^ ^ JztiA) Zf -to Jzu{A) Zi -Uq 



where to = t|Q2=o and Uq = m|q2=o- Thus 



JO 



1 /•! dia dzt {-tf ^ 1 d^b dz^htf 



S Jx+ - Xu Jzt{Al) zf -to S Jx^ ^ - Xt J zt{A2) zf -to 

The evaluation of "+" functions in this basis is quite straightforward. Converting the 
phase space into z^, proceeds along similar lines. 
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